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The monograph discusses the mathematical 
princtples of certain alloved variants of iner- 
ttal guidance of the flight of ballistic missiles, 
t.e., the control without the use of any external 
tnformation (radio signals, radiation of stars 
and others). 


It is assumed that the providing of the 
assigned flight range of the missile ie produced 
as a result of the well-timed ewttching off of 
ite engine by a signal entering from the computer. 
Fed to the input of the latter are current read- 
inge of sensitive elements of the system of iner- 
tial guidance, which measure the apparent accel- 
eration of the missile or integrale from the 
apparent time acceleration. ‘he control should 
be such that the deviation of the actual motton 
‘of the missile from the rated does not have an 
effect on the range of ite flight. For this mea- 
suring instruments of the system of inertial 
navigation are placed aboard the rocket in a 
quite definite manner, and, specifically, the 
direction of their axes of sensitivity tn a num- 
ber of cases is stabilised by meane of gyroscopes. 
The setting of the necessary directions of axes 
of sensitivity and the search of relations which 
should satisfy current readings of sensing ele- 
ments for the formation in the computer of the 
command for the switching off of the engine, and 
this is the main content of the book. Questions 
of the inertial guidance of motton of the missile 
in a lateral direction are also constdered. 


The monograph te intended for specialists in 
the field of the theory of control processes. It 
can prove to be useful in the investigation of 
new probleme of this discipline and also as a 
training manual. 
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THE AUTHOR'S COMMENT 


This book came about as a result of thoughts of many years cover 
the general problems of inertial guidance systems of moving objects. 
Some of them are presented in my other monograph, published in 1963 
under the name of "Mechanics of Gyroscopic Systems," in which, in 
accordance with its name, considerable space has been given also to 
the theory of sensing elements of systems of inertial guidance — a 
‘gyroscope in a gimbal suspension, gyroscopic stabilizer and gyroscuptc 
integrator of the apparent unit accelerations. In this book, on the 
contrary, studied basically is the possibility of the solution of 
problems of inertial guidance of ballistic missiles by means of 
different combinations of gyroscopic and other mechanized sensing 
elements and computers. Of course, questions examined here exhaust 
by no means the whole diversity of ideas and possibilities of auto- 
nomous control of ballistic missiles. In recent years in the theory 
of specific inertial systems great successes have been achieved. 
They have not received as yet the proper reflection in monographic 
literature. It can be hoped that for the assimilation and develop- 
ment of new ccencepts of inertial control the proposed book proves to 


be useful. 


The author thanks A. S. Kachanov, D. M. Klimov, D. F. Klim, 
I. S. Kovner and M. Ye. Temchenko for their valuable advice in the 


editing and preparaticn of the manuscript for press. 
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INTRODUCTION 


Given in this small monograph are mathematical bases of one of 
the possible aspects of the theory of flight control of ballistic 
missiles without any external information (radio signals, radlation 
of the stars and others). The initial data for such control are 
readings of instruments located aboard the missile. They record 
both the crientation of the missile with respect to fixed directions 
(at fixed stars) and the difference in its motion from free fiight 


' without resistance of the medium. The effect of such instruments 
is based on the phenomenon of inertia. Because of this the mentione 
_ guidance is called inertial. Besides instruments which incorporate 


laws of mechanics, in systems of inertial guidance the application 
of instruments founded upon other physical principles (spins of 
elementary particles, standing waves of coherent radiation of lasers 
and others) is possible. | 


Creation of the variant examined in the monograph of the theory 
of inertial guidance is based on the assumption that aboard the 
missile there is stabilization (specifically, gyroscopic) of direc- 
tions of axes of sensitivity of meters of apparent accelerations — 
newtonmeters or instruments directly measuring the time integral 
from current values of apparent acceleration. The latter are called 
below integrators of accelerations; they can also be named impulso- 
meters [Translator's note: This term is not verified]. The selec- 
tion of these sensing elements was determined by the many years of 
interest of the author on gyroscopic instruments and other mechanica 


measuring devices. 


FIDAMT-24-291-79 viii 





i 














The instruments have a different kind of error. The effect of 
the latter on the accuracy of flight of the missile requires a 
Special investigation, which gces beyond the limits of this monograph. 
Here it is assumed that instruments of systems of inertial guidance 
are ideal, i.e., they function as though without errors. The same 
is referred to servodrives and reading devices’ and also any, which 
accomplish, specifically, arithmetical operations and the integration 
of current values of different parameters. 


Airborne instruments and devices should not be excessively 
complex and bulky. Therefore, keeping in mind the use of simplified 
schemes, with the solution of problems of inertial guidance certain 
fundamental errors can be allowed. The latter, however, should not 
involve great deviations in the missile from its target. The selec- 
tion of an appropriate compromise (at which deviations are small, 
and the system of inertial control is not too complex) constitutes 
one of the basic tasks of the designer of a specific system of 
inertial guidance. 


Independently of the form of the guidance system, in its design 
one should proceed from the fact that the same target — the end of 
the free-flight section — can be hit by the ballistic missile, with- 
out necessarily moving along the programmed trajectory but according 
to an infinite set of other adjacent trajectories. Therefore, for 
an accurate hitting of the target it is not necessary that at the 
instant of termination of the active flight section of the missile 
its coordinates and component velocities in a certain system of 
coordinates, fixed relative to the earth, would be equal to the 
calculated provided by the program of control. The last observation 
is very important since the providing of the just mentioned equalities 
places before the system of control of flight of the missile very 
difficult, not always feasible tasks. The attaining, a rather 
accurate determination of current values of the coordinates them- 
selves of the moving missile and projections of its velocity is 
considerably simpler. It is possible, therefore, to interrupt the 
powered-flight section of the missile at exactly that instant when 
the totality of deviations in its moving coordinates and projections 
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of velocity from the appropriate calculation values provides sub- 
sequent motion along one of the trajectories leading to the target. 
The determination of relaticnships which should be satisfied by the 
mentioned deviations at instant of the switching off of the engine, 
constitutes the main task of the stated theory of inertial guidance 
of flight of ballistic missiles. In this case without apparent 
deterioration of the accuracy of hitting of the missile on the target, 
considerable simplifications of guidance systems can be produced 
because of the separate control of the range and lateral motion of 
.the missile and also due to the selection of specific directions of 
axes of sensitivity of inertial instruments and the forming of 


signals from them. 


With gyroscopic stabilization of axes of sensitivity of newton- 
meters and integrators of accelerations, it is natural to construct 
caiculated equations and equations cf inertial guidance in the 
system of coordinates not taking part in the rotation of the earth. 
Let us note, nowever, that range of flight of the missile is deter- 
mined in this case not only by values of its coordinates and com- 
ponents of velocities relative to such a system at the instant of 
termination of the power-flight section, but also of duration of the 


latter. 


* 


In the computer of the given system of inertial guidance of. 
flight range according to current readings of newtonmeters and 
integrators of accelerations, there is generated a certain alternating 
magnitude, uniquely connected at each instant with the magnitude 
of range of the missile. Understood by the latter is that magnitude 
of range which is obtained if the thrust of the engine at the instant 
of time suddenly becomes zero. To provide accurate hitting of the 
missile at the target, the switching off of the engine should be 
produced at the instant of time of achievement by the mentioned 
magnitude of assigned value. Since this magnitude is determined by 
the actual course of the change in coordinates and projections of 
the velocity of the missile, then it is called functional. As it : 
will te shown below, it is possible to propose different forms of 
functionals. With their construction in the computer, subsequently, 
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there will not be taken into account the squares and prcducts of 


deviations in the actual coordinates and projections of the velocity 
of the missile at the instant of end of operation of the engine from 
the rated. The same refers to deviations of moving coordinates and 
projections of the velocity of the rocket and also to the time 
interval between the necessary and calculated instants of the switch- 
ing off of the engine. The assignment of the control of range is 
thus solved in a linear approximation. This should not lead to 
great errors in flight at the assigned range of the missile with 
well-controlled thrust of the engine and with the possibility of 

the switching off of the latter after the feed of the appropriate 
command with a minimum delayed pulse. In the case when control of 
the thrust is hampered, subsequent refinements of equations deter- 
mining the termination of the power-flight section are necessary. 
These problems and also the problem of inertial control without the 
switching off of the engine are not considered in this monograph. 


Results of investigations, given in this monograph indicate 
that in the flight of a missile over great ranges it is expedient 
to have on the gyroscopic stabilizer two integratcrs of acceleration. 
Axes of sensitivity of these integrators shoula ..« definitely oriented 
in twe different directions. Fed into the computer, which determines 
the instant of the switching off of the engine, are the current value 
of one integrator and the result of complementary integration of 
readings of the other. In principle, it is possible to bypass one 
integrator of acceleration if in the process ef flight of the missile 
properly both parameters of the integrator and direction of its axis 
of sensitivity are changed. 


The use of two stabilized integrators of accelerations or one 
but with alternating parameters allows constructing the so-called 
complete functional most accurately resolving the problem of the 
control of range of the missile. Subsequent development of the 
sheory leads, specifically, to the basis of the arrangement of 
integrators of accelerations directly aboard the missile under the 
condition of introduction inte the guidance system of the so-called 
standard integrator of accelerations with accurate execution of its 
commands. 

xi 
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The acceleration of force of the earth's gravity depends on the 
distance between the missile and the center of the earth and also 
on geographical coordinates of the missile. Consequently, forces of 
gravity acting on the missile in its actual and calculated motions 
are net equal. This fact can condition the error in the determining 
of the necessary instant of switching off of the engine even during 
accurate operation of the integrators of accelerations and of gyro- 
scopes of the system of inertial control. For the complete elimina- 
tion of this error the airborne computer should be supplemented by 
a certain comparatively complex integrating equipment. However, a 


_considerable decrease in this error can be achieved also by simpler 


means — because of a small change in the orientation of the integrator 


of acceleration. 


Deviation of the rocket from the target in a lateral direction 
can be eliminated with the help of introduction into the system of 
inertial guidance of the so-called lateral integrator, which controls 
the motion of the rocket in the direction of the normal to assigned 
programmed plane of its flight. The distinction of the actual 
magnitude of duration of the free-flight section from its calculated 
value leads to additional lateral deviations of the missile because 
of the rotation of the earth. However, in principle it is possible 
to avoid errors of such a kind if we use current values of integrators 
of acceleration of system of range control. 


Content of book is the following. Chapter I gives a solution 
to the problem of inertial guidance of ballistic missiles in a 
simplified formulation for the purpose of explaining basic ideas and 
methods examined in detail in subsequent chapters. Here the earth 
js taken as being flat and not rotating. The effect of the atmosphere 
is not considered and the force of gravity is considered constant in 
magnitude and direction. It is natural that under such simplifying 
assumptions the motion of the missile in the free-flight section of 
its flight is completely described by known equations of theoretical 
mechanics about the motion of a material particle in a vacuum in a 


i 
} 
i 
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uniform field of the force of gravity. Further it is explained 
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which relationship should be satisfied by small changes in parameters 
of the end of the power-flight section, i.e., coordinates of the 
missile and projections of its velocity at the instant of the switch- 
ing off of the engine in order that it hits the assigned target with 
an error of not more than the second order of smallness. This 
relationship serves as the basic point for the formatior of a certain 
variable, called the ballistic function. In the examined case the 
value of this function at any fixed instant of time is the approximate 
expression for the error in the range of the rocket, which will arise 
if the engine of the missile is turned off at precisely this instant. 
It is natural that for the hitting of the missile on the target the 
engine should be turned off at the instant of the passage of current 
values of this function through zero. 


To construct a ballistic function aboard the rocket in the form 
of a certain electrical or mechanical magnitude with the variant of 
the inertial system selected in the monograph the presence of special 


“sensing elements is proposed ~ newtonmeter and a computer, whict. 


contains in its composition integrating, multiplying and summing 
elements. In the example of the simplest meter of apparent accelera- 
tion ls established the connection between the magnitude being mea- 
sured, the real acceleration of the missile and acceleration of the 
force of gravity. Concepts of apparent velocity and apparent .travel 
of the missile are introduced and it is indicated, specifically, that 
for the direct measurement of the apparent velocity it is necessary 
that the axis of sensitivity of the integrator of accelerations be 
stabilized. It turns out that the ballistic function, which is 
referred to the given case of motion of the missile in a uniform 
field of gravity without allowing for resistance on the side of the 
atmosphere, can be constructed aboard the missile, using only the 
integrator of acceleration with subsequent integration of its current 
readings by means of the computer. At the end of the chapter an 
analysis of this question is given both from an analytical and 
geometrical point of view. It is indicated that the flight range of 
the missile will be changed by a magnitude of the second order of 
smallness if the vector of its actual velocity at the end of the 
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power-flight section appears the same as that in the calculated case, 
and the position of the missile will be somewhat displaced in a cer-~ 
tain definite direction, or, on the contrary, with a fixed position 
of the missile the difference between the real and calculated vectors 
of its velocity at the same instant of time proves to be definitely 
directed. In the examined case the mentioned directions coincide 
with each other. They, specifically, are parallel to the vector of 
velocity of the missile at instant of its hitting of the target with 
calculated motion. In general when the curvilinearity of the form 

of the earth, the effect of the atmosphere on the descending branch 
of the free-flight section of the trajectory and heterogeneity of 

the field of the earth's gravity are taken into account, such 
directions also exist, but they are not parallel. The perpendiculars 
to them, the so-called \- and u-directions can serve, as it is " 
indicated in more detail in Chapter III, for orientation of axes of 
sensitivity of integrators of accelerations of the appropriate system 
of autonomous control of the flight range of the missile. 


Chapter II gives ballistic functions of various kind for the 
general case of flight of the missile, i.e., without additional 
simplifying assumptions, which took place in Chapter I. Here it is 
impossible to express the magnitude of range of the missile in the 
form of a definite equation, which would contain parameters of the 
end of the power-flight section, for example in the starting system 
of coordinates. Nevertheless, correct to smallness of the .second 
order relative to differences between current parameters of the real 
motion of the missile and parameters of the end of the power-flight 
section of its calculated motion, it 1s possible (with the help of 
the use of the same instruments as in Chapter I) to determine the 
error in the range which would take place with the switching off of 
the engine at the current instant of time. The expression for such 
an arbitrary error in range can be accepted as the ballistic funce 
tion. Derivatives of the magnitude of flight range of the missile 
entering into this function, according to parameters of the end of 
the power-flight section, i.e., according to coordinates and pro- 
jections of velocity, are taken at calculated values of the latter. 
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Thus they are constant quantities, which can be determined earlier 
for each assigned case of flight of the missile. 


With the inertial control of range the current value of the 
ballistic function should be determined onboard the missile itself. 
This leads to the necessity of continuous determination by means of 
the computers of current coordinates of the missile and projections 
of its velocity according to readings of newtonmeters or, in other 
cases, integrators of apparent accelerations. The starting system 
for this goal proves to be barely adequate because of the necessity 
of the calculation of translational and coriolis accelerations. 
Therefore, it is more expedient to pass from coordinates and pro- 
jections of velocity of the missile in the starting system of 
coordinates to appropriate magnitudes referring to the nonrotational 
system with the origin at the center of the earth. As a result, 
after the rejection of terms of the second order of smallness and 
terms dependent on the motion of the missile in a direction per- I 
pendicular to the plane of flight, the so-called initial ballistic : 
function is formed, and they are linear with respect to the moving : 
coordinates and projections of velocity of the missile and time. The 


.-expression for the initial ballistic function in a nonrotating system 


of coordinates contains a derivative of the magnitude of flight range 
of the rocket with respect to duration of the power-flight section. ] 
However, this derivative can be excluded from the analysis by means 
of the use of some auxiliary relation. The latter follows from the 
equality of the actual range to its calculated value, if coordinates I 
and projections of velocity of the end of the power-flight section 
accuratély coincide with current values of coordinates and projections 
of velocity of any particles of the calculated free-flight section. 
Finally the basic ballistic equation is obtained for deteorr*rotion 

of the necessary current time of the switching off of the engine in 
flight of the missile at the assigned range. Determination of 

current coordinatés and projections of velocity of the missiles, 

which enter the left part of this equation, requires continuous 
integration aboard the missile of appropriate nonlinear differential 
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equations, which connect the second derivatives of coordinates with 


"a 
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appropriate projections of the apparent acceleration and accelera- 
tion of the force of gravity. However, we can substantially simplify 
such a problem, if into the basic ballistic equation we substitute 
expressions of coordinates and projection velocities obtained as a 
result of the integration of approximate linear differential equations 
in the so-called isochronal variations of coordinates of the missile. 
The latter are differences between the acutal and calculated values 
of appropriate coordinates of the missile, which refer to the same 
instant of time. With’ the proper arbitrary extension of the cal- 
culated ‘power-flight section of flight of the missile during the 
calculated instant of the switching off of the engine isochronal 
variations of coordinates, projections of velocity and projections 

of apparent acceleration can be considered small magnitudes during 
the whcle interval of time which corresponds to the real power-flight 


. section. 


Differential equations. for isochronal variations of coordinates 
are derived in Thapter III, and their approximate integration, allow- 
ing ror terms appearing as a result of the heterogeneity of the field 
of the earth's gravity, is developed in Chapter IV. In Chapter III 
in differential equations for isochronal variations of coordinates, 
the mentioned terms are completely droped, whereupon variations of 
projections of velocity become equal to variations of projections 
of the apparent velocity of the missile, and variations of the 
coordinates themselves — to variations of projections of the apparent 
path. This allows converting the basic ballistic equation to one of 
the forms allowing the construction of current values of its left side 
directly onboard the missile according to current readings of integra- 


tors of accelerations. 


Different orientation of axes of sensitivity of the integrators 
of accelerations are possible. Depending on this, the computing . 
part of the system of inertial control of the range and the number 
of elements entering it are changed. For example, with orientation 
of axes of sensitivity according to invariable \- and u-directions, 
which was already mentioned above, it is possible to manage with only 
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one element of repeated integration of accelerations unlike the case 
of using two such elements necessary in the arrangement of axes of 
sensitivity in parallel to the axes of the nonrotating system of 
coordinates. 


It proves to be possible to construct a ballistic equation even 
with the help of a single meter of acceleration, if we according to 
the earlier assigned law change the orientation of its axis. In the 
computer there should be provided in such a case the integration of 
current readings or this meter, preliminarily multiplied by the 
coefficient, whic1 also is changed with the course of time. 


The problem indicated above on the creation of a ballistic 
equation can be solved approximately also with the help of one 
integrator of acceleration, as this is indicated at the end of the 
chapter in the example of construction of the system of inertial 
control of range with accurate control of the direction of thrust 
of the engine according to indications of standard integrator of 


accelerations. In this case the axes of sensitivity of the standard 


integrator and integrator of the system of control of the range are 
located perpendicular to each other. Finally, if we consider the 


. deviation in the force of thrust of the engine from the longitudinal 


axis of the missile to be insignificant, then for solution of the 
problem of control of range one can use current readings of the 
integrator, the housing of which 1s directly fastened to the side of 
the missile, and the axis of sensitivity is parallel to its longi- 
tudinal axis (so-called longitudinal integrator). 


In Chapter IV there is formed the integration of the totality 
of differential equations for isochronal variations of coordinates 
of the missile with the approximate calculation in them of terms 
which characterize changes in projections of the force of gravity 
during motion of the missile according to the law different from 
the calculated. Corresponding to this, the basic ballistic equation 
derived out in Chapter II is converted to the form similar to that 
examined in Chapter III. Because of a specific method of integration 
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of the mentioned totality of equations, the problem on inertial con- 
trol of the range is solved by the same means as it was earlier, i.e., 
with the help of current readings of integrators of accelerations and 
computers, which contain only integrating, multiplying and summing 
elements. This method is based on the preliminary replacement of 
equations of the mentioned system by equivalent integral differential 
and cumulative relations with the subsequent introduction into the 
terms being integrated instead of the sought functions, which repre- 
sent current values of coordinates of the missiles approximating 
their polynomials of the second or third degree. Coefficients of 
approximating polynomials of the second degree are selected so that 
the polynomials would satisfy the same initial conditions as the 
sought functions, and, furthermore, coincide with the latter for the 
instant of termination of the power-flight section. With a more 
accurate solution to the problem, one should approximate the isochronal 
variations of coordinates by means of polynomials of the third degree. 
‘In this'case time derivatives are equated also to appropriate pro- 
jections of velocity of the missile at the instant of termination of 
the power-flight section. As a result for the determination of 
variations of coordinates and projections of velocity of the missile, 
algebraic equations with coefficients dependent on time are obtained. 
The latter, with sufficient accuracy, can be replaced by constants 
equal to values of these variations at the instant of termination of 
the power-flight section of flight of the missile. Finally the 
desired variations are expressed in terms of readings of integrators 
of accelerations and integrals of their time readings. Thereby, the 
subsequent creation of different forms of ballistic functions similar 
to those given in Chapter III, fundamental has no basic difficulties 
and is reduced basically to certain changes in coefficients and 
directions of axes of sensitivity of integrators of accelerations. 


Finally, Chapter V of the book is devoted to the general 
problems of inertial control of the lateral moticn of ballistic 
missiles. The removal of lateral deviation in the rocket from the 
target requires different methods than the providing of the ascigned 
range of its flight, where the required accuracy can be obtained 
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because of the well-timed switching off of the engine. The duration 
of the power-flight section with flight at the same range depends for 
rockets of the same design on a number of random facts and, con- 
sequently, is itself a random magnitude, the mean value of which is 
near to the magnitude of duration of the power-flight section of the 
calculated motion. Therefore it is sufficient that the control of 
lateral motion of the missile would be especially accurate only during 
the interval of time in which switching off of the engine can occur. 


The current removal of the missile from the programmed plane is 
connected with projections on the normal to this plane of its apparent 
acceleration and acceleration of the force of gravity by means of a. 
differential equation of the second order. To control the lateral 
motion of the missile, it is necessary to know current values of 
this removal and its time derivative. They can be obtained by means 
of the so-called lateral integrator of accelerations and simplest 
computer, As a first approximation can in the mentioned differential 
equation drop terms containing the projection of the force cf 
gravity, whereupon it is immediately integrated. The lateral depar- 
ture of the rocket proves to be equal to the projection of its 
apparent path on the normal to the programmed plane of flight, and 
the value of velocity of this departure — respectively, to the pro- 
jection of the apparent velocity. With the obtaining of more accurate 
expressions for the mentioned sought magnitudes, the approximation 
method of integration of differential equations given in Chapter IV 


is used. 


In Chapter V it is also indicated how on: can use readings of 
the instruments of inertial control of the renge in order toc remove 
the additional lateral deviation in the rocket from the target because 
of noncuincidence of its actual and calculated motions in the pro- 
jection on the programmed plane. 


In the appendix of the book an analytical derivation of cosines 
of angles between axes of system of coordinates with a finite turn 
is given. The table of cosines is used in Chapter II. 
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CHAP? ER I 


THEORY OF INERTIAL CONTROL OF FLIGHT RANGE OF THE 
BALLISTIC MISSILE IN SIMPLIFIED FORMULATION 


§ 1. Equation Expressing Flight Range of the Missile 
in Terms of the Coordinate an rojecticn of 
Its Velocity at the End of e 


ower-Flight Section 


The general formulation of the problem of control of Tlight 
range of the ballistic missile is most convenient to exp sin iv the 
following simplest example. 


Let us assume that the earth is flat and does not rotate, the 
atmosphere surrounding the earth is absent,*and the acceleraticr cof 
the force of the earth's gravity f, coinciding in this case with 
gravity g, is constant in magnitude and direction. Equations of the 
motion of the missile’ on the free-flight section of its trajectory 
under such simplifying assumptions are completely integrated, and 
thereby the whole solution to the problem can be conducted in general 
form up to the end. 


Let the power-flight section of flight of the missile proceed 
from the origin 0 of the fixed system of coordinates zy, the x axis 
of which is horizontal, and the y axis is directed along the vertical 
upward (Fig. 1). The initial veiocity of the missile is considered 
equal to zero. 


lMHere and further the term “motion of the missile" is understocd 
as motion of its center of mass. 
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Fig. l. 


The caleulated or programmed motigqn of the missiles is that one 
at which its coordinates are changed in accurate conformity with their 
precalculated values. 


Let us designate coordinates of the real motion of the missile 
at instant o, termination of the power~flight section of its flight, 
i.e., at. the instant of complete switching off of the engine, by zx 
and y, and the same coordinates in calculated motion — by z* and y* 
(Fig. 1). Let us designate by v, and Yy projections on the axis x 
and y of the velocity of the rocket at the instant of the termination 
of the power-flight section of its real motion and respectively by 
ve and v* — their calculated values. The magnitude of duration of 
the pcower-flight section of the actual flight of the missile o in 
generai is distinguished from its calculated value o®, 


The integration of equations of motion of the missile in the 
pow... Flight section of motion 


SPao, Bus, (1.2.2) 


immediately leads to the following expressions for current projec- 
tions of velocity v(t) and Bers 
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v(t) = 20 ov, vy (t) = “3 oo, — gt, (1 1 2) 





and its moving coordinates x(t) and y(t): 
: < _ 
z(t) mz + og, y(t) yt of——-. (510 35 


Time t is read off here anew, i.e., from the instant of termination 
of the power-flight section and, consequently, from the beginning of 
the power-flight section of motion of the missile. Arbitrary con- 
stants with integration are selected so that there would take place 


the apparent equalities 
vo, (0) =O, Py (0) =p = (0) az, y (0) my. ( 1 : 1 7 4 ) 


Considering that the free-flight section is finished at instant 
t = tT at point C ~— intersection of the trajectory of the missile 
with the x axis, it is easy to determine the connection between the 
range of its flight Z and appropriate parameters zx, y, v, and vp. — 
end of the power-flight section. Setting, first of all, in the 
second equation (1.1.3) t¢ = t and y(t) = 0, we obtain the following 
quadratic equation for determination of the duration of the free- 


flight section Tt: 


get — 20,8 ~ 2y =. (1.1.5) 
Hence . 
tm (0, + Voy + 289), (1.1.6) 


since one takes the positive root of equation (1.1.5). Having 
substituted now the found value t into the first equaticn (1.1.3), 
we obtain the desired expression for the range of the rocket as 
functions of parameters zx, y, vS and vy -— end of the power-flight 
section, namely: 


bm 5(t) = 24-2 (0, + Voy + 265). Cleke7) 
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§ 2. Condition of Invariability of Flight Range with 
Little Distinction in Coordinates and Projections 
of Velocity of the Missile from Their 
Calculated Values at the Instant 
of Termination of the 


Power-Flight ° 


Section 


Let us assume that the actual values of parameters of the power- 
flight section of flight of the missile zx, y, », and y are distin- 
guished from their calculated values by small magnitude 


Acmr—z, Ayay—y, Av, = 0,— 05, A, = ¥,— vy. (1.2.1) 


The flight range of the rocket Z, in accordance with that given in 
§ 1, and, specifically, according to equation (1.1.7), is a function 


of magnitudes zr, y, v, and Py» i.e., 


I = E(x, y, Var 0). (1.2.2) 
In accordance with equality (1.2.1), we have 
Tmt (z Az, y+ Ay, 02 + doy 0, + Ady). (49233) 


Expanding now the right side of equation (1.2.3) into Taylor series 
for the function of many variables, we obtain correct to smallness 
of the first order relative to magnitude Az, Ay, Av, and AQ 


eo eee at a of a 
Pm l(z yen or + a As + 5 Ay +57 Ave + 5, Any (1.2.4) 
Partial derivatives of function l(zx, y, Dog ) with respect to 


varlables x, y, Lae and %y are themselves functions of the same 
variables. In this case they can be presented in evident form as 
aresult of the differentiation of function (1.1.7) according to 
he appropriate variable, namely: 


a at. y at 4 3S, 
a Vaan H=-+(,4+V item). 


al 9 ’ ; 
mm 2/44 ——%— }. 
ee, 7( +7z4=) (1.2.5) 
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According to rules of the composition of Taylor series, in 
equation (1.2.4) arguments of aforementioned partial derivatives 
should be assumed correspondingly equal to magnitudes zr*, y*, vi 
and v*, Thus they are for the selected calculated motion of the 
missile certain constants. Taking into account, furthermore, that 


expression 


I ww l(" ¥'. Ben 95) : (1.2.6) 


is the calculation flight range of the missile, let us transform 
equation (1.2.4) co the form 


% tye : 
SS iet yea YS t wt Vy + 2ey ) av.+ 


where the difference 


is the change in flight range of the missile because of noncoincidence 
of magnitudes of parameters of the end of the power-flight section 
Ly Yo VU, and vy with their calculated values zr‘, y*, ve and ie 

The expression for Al, allowing for equation (1.1.6), can be 
represented now in the form 


Al = As + kAy + 1°(Av, + kAD,), (1.2.9) 


where 
e 
%e 


"Trea 


(1.2.10) 


Let us note, that coefficient k has a simple geometric meaning. It 
is easy to show, using equations (].1.2) and (1.1.€), that it is 
equal to the magnitude of the angular coefficient of the normal v to 
the calculated trajectory of motion of the missile at point Cc, 

its intersection with the x axis (Pig. 1). Thus, 


nce 


vO 





bautga, (1.2.22) 


where a — angle between the mentioned direction v and = axis (and 
also between the vertical line and vector of velocity of the missile 
at the instant of its incidence on the earth). 


From equation (1.2.9) it follows that correct to smallness of 
the first crder inclusively the range of the missile remains the 
very same if changes in parameters dz, dy, dv, and Av. — end of its 
power-flight section in comparison with their calculated values, will 
be subordinate to the condition 





An + bay + ¥° (An, + kAp,) = 0. Cio 125 


It is important that this condition is the linear relation with con- i 
stant coefficients known earlier for the selected flight of the : 
missile. 


The duration of the power-flight section of fiilght o, in con- 
trast to parameters cf its end z, y, vs and v_ , does not play any 
role in the determination of the magnitude of range 7. Therefore, 
there is no meaning i> strive with execution of the specific program 
of motion of the missile for the rated value of duration o* of the 
power-flight section. On the contrary, by slightly lengthening or 
shortening the power-flight section, it is possible to obtain a 
fulfilling cf just the given condition (1.2.12) and, consequently, 
provide the calculated range of the missile at not exactly an accurate 
sequence of the actual coordinates to their programmed values. 


§ 3. Error in Flight Range of the Missile Presented 
in the Form of a Function of the Duration 
of the Power-Flight Section 


Let us assume that, as also in § 1, z(t) and y(t) are functions 
which are the actual (not calculated) change in time of coordinates 
of the missile on the power-flight section, a vi(t) and ult ts re- 
spectively, are projections of its velocity, and the instant t # 0 
at this time corresponds te the beginning of the power-fiignt section. 


€ 





Let us substitute into equation (1.2.9) for Al quantities Az, 
Ay, Av, and Av, respectively by differences 


5), WO—¥, Oa) —02, 0, (— 04. (1.3.2) 


As a result let us obtain the time function 
a(t) w2()—2 +k ly(t)—y'] +o [ve (t)— 92] + 
+ ke’ [vy(t)— 051, (15362) 


which we call baltlisttie function. 


At'the instant of the switching off of the engine t = of this 
function, according to equation (1.2.9), turns into AZ and, con- 
sequently, determines the error in flight range of the rocket because 
of an accurate execution of the program of the power-flight section 
but accurate observance of its duration. However, if we interrupt 
the power-flight section at that instant when the dallistic function 
e(t) turns into zero, then the error into the range proves to be a 
magnitude of the second order of smallness, and the desired accuracy 
of flight of the missile will be achieved. Hence it follows that 
the problem of inertial control of the flight range of the missile 
can be reduced to the construction of current values of the ballistic 
function directly aboard the missile and to the switching off of its 
engine with the passage of function e(t) through zero. For this 
purpose, besides the summing and multiplying devices, the presence 
aboard the missile of special instruments which measure the apparent 
acceleration with subsequent integration of their current readings, 
or instruments directly recording the time integral from the apparent 
accelerations — integrators of accelerations is necessary. Meters 
of apparent accelerations will subsequently be called ewtonmeters. 


§ 4, Connection Between the Acceleration of the Missile 


in Its Motion Relative to the Eixed System of 
Coordinates and Reading o 


Newtonmeter Set on It 


In the system of control of the flight range of the ballistic 


o 
i 


missile when using newtonmeter one should keep in mind that the latter 
can measure only the apparent and not the real acceleration of that 
place of the missile where they are located. The apparent accelera- 
tion is usually called the difference petween the acceleration of 

any points in the fixed system of coordinates and the acceleration 

of the force of gravity. Specifically, standard uniaxial newton- 
meters should measure the projection a, of the apparent acceleration 
on its axis of sensitivity v, 1.e., magnitude 


Game, (1.4.1) 


where We projection on axis v of the real acceleration of the 
newtonmeter relative to the fixed system of coordinates and I ~ 
projection on the same axis of acceleration of the force of the 


earth's gravity. 


Let us explain equation (1.4.1) in the example of the newton- 
meter, tre sensing element of which is a small weight of mass m 
attached to the end of the spring of rigidity ec (Fig. 2). The other 
end of the spring is sealed in the housing of the newtonmeter. The 
smali weight can be moved within the housing on straight line v, 
which is the axis of sensitivity of such a meter of apparent accelera- 


tion. 





Fig. 2. 











The equation of motion of the small weight relative to the 
housing of the newtonmeter, if we disregard the mass of the spring 
and friction cf small weight about the internal cylindrical surface 
of the housing, can be presented in the form 


~ 


—m Se t+ mul a mf, +d. (1.4.2) 


« - 


Here §6 — movement of the small weight from that position E-«at which 
the spring has not been stretched; we > projection on direction v cf 
tranctational acceleration, i.e., acceleration with respect to the 
fixed system of coordinates of that place cf the housing where at the 
giver instant the small weight is located (projection of the coriolis 
acceleration to direction v is equal to zero). 


With translational movement of the housing of the rewtonmeter 


tow. Chel 3) 


“where w — acceleration of the piace of fastening of the spring to 
the housing. 


If, however, the housing of the newtonmeter has, furthermore, 
angular motions, then equality (1.4.3) should be considered as 
approximate. However, because of the comparatively small dimensions 


of newtonmeters the difference between we and wy is important. 


Let us assume that the frequency p of natural oscillaticns of 
the small weight, determined cy equation 


p= V<. (1.6.4; 





'It 4s possible to show that this difference is equal tc ‘he 
product wee where wy — projection cf angular velocity of ‘the nousing 
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is sufficiently great. Then the amplitude of these oscillations, 
conditioned basically by the change in acceleration of the housing 
of. the newtonmeter w. will be small. In this case the term -ms 

in equation (1.4.2) can be dropped, and allowing for equality (1.4.3) 


we can obtain relation 
man, ax mf, + 8. (1.4.5) 


The member cé of equation (1.4.5) is the elastic force of the 
spring. Deformation 6 (see Fig. 2) will be considered positive if 
the small weight is displaced in a negative direction of axis v. 
Wher 6 = 0 the spring is not stretched, and the small weight is 
located at position F& (see Fig. 2). The observable shift in the 
small weight from position EF can be graduated so that it would 
directly measure a certain magnitude a, connected to the deformation 


of the spring by relation 


= 38. (2.4.6) 


Replacing here the magnitude of deformation of the spring 6 by its 
expression, by tne following from equation (1.4.5), we will arrive 
at equation (1.4.1). Thus, the shift in the small weight is pro- 


portional to the projection of its apparent acceleration on the 
direction v, 1.e., on the axis of sensitivity of the newtonmeter. 


The current reading of the integrator of the apparent accelera- 
tions or simply integrator of accelerations (it can also be called 
tmpulsometer) is a time integral from the projection of the apparent 
acceleration on its axis of sensitivity, 1.e., magnitude 


t 
Vo(t) = feed. (1.4.7) 


If the axis of sensitivity of the integrator of acceleration 
retains a fixed direction with respect to the fixed system of 


coordinates, then 
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Oe, (1.4.8) 


where vft) — projection on the direction v of velocity of the hous- 
ing of the integrator in the same fixed system. Therefore, sub- 
stituting into equation (1.4.7) the expression for a.) from equation 
(1.4.1) and using equality (1.4.8), we obtain the relation 


? . . 
Vo(t) = 24 (8)— 9, (0) — { fournat, (1.4.9) 


which connects the projection of velocity v(t) with the reading of 
the integrator of acceleration Vij(t). The latter is called in this 
case the projection of the apparent velocity on direction v. 


Function fi(t) which stands under the sign of the integral in 
the right side of relation (1.4.9) when the earth is not proposed to 
be flat, is changed with the course of time because of a change in 
the position of the missile relative to the earth. 


With alternating crientation of the axis of sensitivity v, for 
example, when the integrator of acceleration is located directly 
aboard the missile, equations (1.4.8) and (1.4.9), of course, are 
not applicable, and the reading of the integrator of acceleration is 
no longer equivalent to the corresponding projection of the apparent 
velocity. 


§ 5. Construction of a Ballistic Function by Means of 
Two Integrators of Accelerations and a Computer 


In this section let us examine the use in the system of inertial 
control of the flight range of two integrators of accelerations. 
Let us assume that the integrators are stabllized so that the axis 
of sensitivity of one of them during the whole power-flight section 
of tne missile would remain parallel to the horizontal x axis, and 
the axis of sensitivity of the other wouid be directed parallel to 


the vertical y axis. 


Ll 
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In the examined case, i.e., under the assumption of a nonrotat- 
ing flat earth, projections of accelerations of the force of gravity 
on the x and y axes are expressed by equations 


fam0,  fym—e. (1.5.2) 


Consequentiy, in accordance with formula (1.4.1), projections of the 
apparent acceleration on the same axis are the quantities 


6, = w;(t), ay mw, (t) +8, : (1.5.2) 
WHEE do, (t) dz (t) dv, (6) d¥y (t) (1 5 3) 
oan, wn y 2, 5. 


are corresponding projections of the real acceleration of the missile 
in the fixed system of coordinates zy. Here x(t) and y(t), just as 
earlier, are current coordinates of the missile on the power-flight 
section of its flight. . 


The current readings of the integrators with axes of sensitivity 
parallel to axes of coordinates x and y, according to equation 
(1.4.7), are quantities 


eesti ine mene oe 


t e 
Van J eattat, Ve = fede, (1.5.4) 


which are projections of the apparent velocity of the missile, 
respectively, on the axis of the fixed orientation x and y. In 
virtue of equations (1.5.2) and (1.5.3) they are connected with 
projections v(t) and v (t) of the real velocity of the missile 
relative to the system of coordinates zy by relation 


Vi(t) = ve(t), Vy (t) oy (t) + at. (1.5.5) 


Into the latter constants of integration are absent, since projec- 
tions of the real veloe¢ity of the missile at the initial instant . 
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of the power-flight section, i.e., when t = 0, are equal to zero. 


Coordinates of the missile x(t) and y(t) at this instant are also 
equal to zero. Therefore, as a result of the integration with 
respect to time of current indications vi) ane vo) of corre- 
sponding integrators of accelerations will be reduced cn the basis 
of equations (1.5.3) and (1.5.5) to the foliowing relation: 


S$10=2(),  H(Namr +H. (1.66 


Tn them time functions 


4 e : 
SAM PVN, — 5, —= SV, (at Cae: 
can be called projections of the apparent path of the missile. 


By means of relations (1.5.5) and (1.5.6) one can express , 
voordinates x(t) and y(t) in terms of functions S(t) and 5, ‘*) 
and projections v(t) and v(t) of its real velocity — in terms of 
current readings of integrators of accelerations Vi(t) and V, (t.. 
If further we substitute these results into equation (1.3.2) for 
the ballistic function e(t), then after the simplest conversions we 
will obtain for it the following expression: 


868) = Salt) + 8S, () + 1V 4 (0) FAY, (91 — 
a byt vs +e) — chet — hE, (1.5.8) 


where the coefficient k, as before, is determined by equation 
(1.2.19). 


For the construction aboard the missile of a current value of 
the ocllistic function eft) in the form of an alternating mechanical 
or electrical magnitude fn the examined case the presence of two 
integrators of accelerations, clccks and computer is necessary. : The 
composition <f the latter, specifically, should include two elements - 
of additionai integration of readings of integrators of accelerations 





§ 6. Control of Fld ht Range of the Missile b Means 
of a Single integrator of Accelerations with 


inclined Orientation of Its Axis 


An attentive examination of expression (1.5.8) for the ballistic 
function e(t) leads to the conclusion that sums 


Va(t) + kV, (2), 5, (t) + 45, (8), (1.6.1) 


entering into its composition can be formed aboard the missile by 
means of unl; the integrator of accelerations with subsequent addi- 
‘tional time integration of its instant readings. The axis of 
sensitivity v cf such an integrator should be inclined during the 
whole time of th2 power-flight section of flight of the missile to 
the horizontal axis x at a constant angle a, the magnitude of which 
is connected witn coefficient k by relation (1.2.11). Thus (see § 2 
of this chapter), the direction of the axis of sensitivity of the. 
mentioned integrator is perpendicular to the vector of velocity of 
‘ne missile at the instant of reaching by it of the target C with 
eniculated moticr. (fig. 1). 


Let us note, first of all (Fig. 3), that for projections on 
the direction v of the acceleration of the force of gravity and real 
acceleration of the missile we have equations 


i. = —gsing, w, = w,(t) cosa + w,(t) sina. (1.6.2) 


Consequently, according to equation (1.4.1), we have 


a, mm w(t) — fo(t) = w(t) cosa + Iw, (t) + g}sina (1.6.3) 


or, taking into account equality (1.5.2), 


@, (t) = @,(t) cosa + @,(t) sina. “1.6.45) 


The last relation indicates that the apparent acceleration is a 
vector whose components are parallel to axes of coordinates x and y, 


respectively equal to quantities a(t) and a(t). 
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Fig. 3. 


Integrating expression (1.6.4) for the projection of the 
apparent acceleration on the direction v with respect to time, and 
taking into account equations (1.4.7) and (1.5.4), we obtain equality 


Vo(t) = V(t) cosa + V,(t)sina, (1.6.5) 


which indicates that quantity v(t) is in this case a projection on 
the direction v of a certain vector with components vift) and ne 
along axes x and y. 


Function v(t) is the current reading of the integrator of 
acceleration whose axis of sensitivity has the direction v. Integrat- 
ing this time function and taking into account equation (1.5.7), we 
obtain a new relation 


S.(t) = S,(t) cosa + S,(¢) sina, (1.6.6) 
in which quantity 


@ 
‘So(t) = Jvourae (1.6.7) 


ean be called the projection of the vector of apparent path of the 


missile on direction v. 


Relations (1.6.5) and (1.6.6) can be given with the help of 
equation (1.2.11) the following form: 
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ie Vi(t) = cosa [V.(t) + AV’, (2)} (1.6.8) 


5S. (t) = cosa [S.(t) + &S, (t)}. (1.6.9) 


From the obtained relations it follows that sums (1.6.1) and, con- 
sequently, and expression (1.5.8) for the ballistic function e(t) 
can be constructed aboard the missile by means of the use of current 


readings of the single integrator of acceleration and of its time 
integral. Really, on the basis of relations (1.6.8) and (1.6.9), 
expression (1.5.&) 1s converted to the form 


0(8) me he 154(0) + V(t) fe" ay") — 
—T (vr + ko) — chgt—k &. (1.6.10) 


As was already mentioned above, the duration of the power-flight 
section is determined by equation e(t) = 0. In accordance with 
equations (1.6.10) and (1.2.11), this equation, which will sub- 
sequentl;, te called ballistie, is reduced to the following, more 
convenient form: 


Volt) +t S.(t) m vecosa + (0, + gt)sina + | 
+4[z* coset (y° +P) sing].- (1.6.22) 


§ 7. Conversion of the Ballistic Equation by Using 
a New Rectangular System of C»ordinates 


The ballistic equation (1.6.11) acquires a more transparent form, 
if we introduce a new fixed system of coordinates En with the same 
origin as that for the system ry. Axis € of this system of coordi- 
nates is directed parallel to straight line v, 1.e., at the same 
angse a to the horizon (to the w axis) at which the axis of sensitivity 
of the integrator of accelerations should be located (Fig. 4). Then 
axis n will be the antiparallel to the vector of velocity of the 
missile at the instant of Llts hitting of the target with calculated 


motion (see Fig. 1). 


bee 
Cts 








Fig. 4. 





Coordinates of the end of the calculated power-flight section 
of the trajectory of the missile in the new system of coordinates 
are designated by £* and n®. They are comnected with csordinates 
x* and y* of the same points in the xy system by equations of con- 


version 


Pas eccsaty'sing, 9a —s*sin's +'y cose. (1.7.1) 


Let us designate projections of velocity of the missile on the axis 
5 and n with its calculated motion at the instant of the switching 
off of the engine by ve and va: They are connected with projections 
of the same velocity on the axes x and y, with quantities us and vie 
‘by similar equations 


my = 9, 0084 + 0, SING, p, = — vs sing + vscoea. (1.7...) 


For moving coordinates of the missile in the system €n, which refers 
to its real motion, let us introduce designations €(t) and n(t-, for 
projections of its velocity on these axes — vp, ft) and vift)s and, 
finally, for projections of acceleration — welt) and w, fe). 


Since the axis £ and direction v are_paraliel, then the oro- 
jection of acceleration of the missile on the direction v is the 


equation 


Ho 
“a 


ww. (8) ow wy (2) = (1.7.2) 
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and, consequently, according to equations (1.4.1) and (1.6.2), we 
obtain 


SP =n +1='e()— going. (1.7.4) 


At the initial instant (t = 0) of the power-flight section, 
coordinates of the missile €/t) and nft) and also projections of 
its velocity v(t) and v, ft) are equal to zero. Therefore, succes- 
sively integrating expression (1.7.4) for d®e(¢) sae? with respect to 
time, we obtain equations 


(ft) = 40 V.() — gi sina, 
&@ = $.()—F sing, (1.7.5) 


where V(t) and S(t), as before, are projections of the apparent 
velocity of the missile and its apparent path on direction v or, 
which is the same, on axis &€. They are connected with the projec- 
tion of the apparent acceleration 2, (t) by relations (1.4.7) and 
(1.6.7). 


When using equations (1.7.5), (1.7.2) and (1.7.1) the ballistic 
equation (1.6.11) for determination of the instant of the switching 
off of the engine is reduced to the following simple fo:m:? 


Hi) + EEE EE. (1.7.6) 


' Not entering into equation (1.7.6) is either the coordinate of 
the missile (t) or projection v, ft? of its velocity on axis n. This 
means that the small cistinetion in the mentioned magnitudes at the 


enw of the power-flight section from their calculated values n*®* and 


'Squation (1.7.6) can be obtained directly from the ballistic 
equation e(t) = 0, in which function eft) 1s taken in the initial 
form (1.3.2). In erder to be convinced of this, it is enough to sub- 
stitute in expression (1.3.2) coefficient k by its value following 
from the relation (1.2.11) and to use equations (1.7.1) and (1.7.2) 
and also those similar to them for quantities &(t) and v(t). 


vi correct to smaliness of the second order should not affect the 
rlight range of the missile. In essence this is explained by the 
fact that for the control of the flight range of the missile it is 
enough to use only the integrator of accelerations whose axis of 
sensitivity 1s parallel to the axis & or, which is the same, as 
direction v. 


§ 8. Geometric Examination of the Condition of 


Hitting of the Missile on the 
Assigned Target 


From a geometric point of view the ballistic equation (1.7.6) of 
the previous section can be given the following interpretation. Let 
us assume that the motion of the missile is such (Fig. 5) that at 
the instant of the switching off of the engine it proves to be at 
point A', shifted with respect to points A — the end of the cal- 
culated trajectory of the power-flight section on a small segment 
parallel of axis n. If its velocity at this instant appears the 
Same as that at the end of the power-flight section cf the calculated 
motion, then the trajectory of the missile in the free-flight section 
will have the same form as that in the calculated case. In order to 
obtain this trajectory, it is enough to shift by the magnitude of 
segment AA’ the whole calculated trajectory forward in the directicn 
of axis n or, which is the same, in the direction of the tangent to 


the trajectory at point C of its intersection with axis x (see Fig. 5). 


It is obvious that the point of intersection of the "shifted" tra- 
Jectory with the x axis will be remote from points Cc, l.e., from the 
calculated points of fall of the rocket, on a smallness of the 
second order. 


Let us assume that, on the contrary, at the end cf the power- 
flight section the rocket arrives accurately at the calculated point 
Ay but at a velocity somewhat distinguichcd from the calculated in 
magnitude and in direction (Fig. 6 and Fig. 7). The form of the 
trajectury proves to be different, and 1: general the rocket will 
nu iunger hit the assigned target. 
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Fig. 7. 


Nevertheiess, as follows from the ballistic equation in the 
form (1.7.6), if the velocity of the missiie will be different from 
the calcuiated by a small vector parallel of axis n (Fig. 8), then 
the deviation of the rocket from the target will also be a magnitude 
of the second order of smallness. Thus, in the case of a nonrotating 


flat earth deprived of an atmosphere, it is possible to show such 


cctneiding directions of smaii vecvors ef the change in velocity or 
the shift if the rocket at the end of the power-flsght section, in 
the presence cf which the change in flight range of the rocket has 
“ho Magnitude of the second oxder of smallness. Similiar dalrections 


be shown in general if one considers the rotation of the earth 
anu its nouspherical state in the presence of an atmosphere. Hewever, 
i 


thev preve to be i fferernt, and for the solution of the problem, on 


26 





selection of the proper instant of the switching off of the engine 

of one integrator of accelerations with the axis of sensitivity of 
constant direction, in general, they are insufficient. The mentioned 
problem (also correct to smallness of the second order) is most 
Simply solved with the help of two integrators of accelerations (see 
§ 3cf Chapter III) with additional irtegration of readings cf <ne 

of them, 
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CHAPTER Il y 


THEORY OF INERTIAL CONTROL OF FLIGHT RANGE OF 
BALLISTIC MISSILES IN THE 
GENERAL FORMULATION 


§ 1. Expression of tne Error in Flight Range of the 
Missile with Small Changes in Parameters of the. 


End of the Power-Flight Section. 
~ Ballistic Function 


The theory sf inertial control becomes incomparably more com- 
olex if with the motion of the missile, in contrast to the simplify- 
ing assumptions accepted in Chapter I, we take account of the change 
in acceleration of the force of gravity both in magnitude and in 
direction, cunsider the earth no longer to be fixed, and in the cal- 
culation of the free-flight section take account of the resistance 
of the atmosphere. The basic difficulty consists here in the selec- 
tion of a certain rather simple function of parameters measured 
aboard the missile by inertial sensing elements. The function 
should be such that with the achievement by it of the earlier 
assigned value, it was possible to produce a switching off of the 
engine, having provided the calculated range of the missile. An 
example of a similar kind was given in Chapter I. This function, 
calle? subroquently, just as in Chapter I, ballistic (sometivcs it 
is called euntrolling and also controlling functional), is con- 
structed aboard the missile by means of the computer, which uses 
current readings of the integrators of accelerations. The magnitude 
of the tallistic function should be directly connected with the error 
in the flight range of the rocket which would oceur if the switching 
off of the engine occurred at the current instant of time. The 
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System of inertial control in turn st:ould give a signal for cessation 
of the operation of the engine upon achievement by the ballistic 
function of a value which corresponds to the turning of the mentioned 
error into zero. 


¥ 


One of the ballistic functions is the errcer itself in deter- 
mining the flight range of the missile, expressed in terms of its 
moving coordinates x(t), y(t), 2(t) and projections vi ft), vers 
vi ft) of velocity relative to the so-called starting system of 
coordinates xyz connected with the rotating earth. In this case it 
is proposed, cf course, that at precisely the instant t complete 
switching off of the engine occurs. 


It is convenient to put the beginning of starting system at the 
center of the earth, having directed the y axis along the radius cf 
the earth through the points of start O and the z axis so that the 
coordinate plane ry, called the program, contains the assigned points 
of fall of the missile ©. If we direct the x axis to the side of 
points ¢, then direction of the z axis is thereby completely deter- 
mined (Fig. 9). It is useful to note, that the trajectory of the 
missile does not lie in plane xy, passing, nevertheless, in the cai- 
culated case through the origin of coordinates U0 and through point °“. 


For the magnitude of the flight range of the missile Zz we can 
take the shortest geographical distance between the point of ‘ 
and actu#l point of fall of the missile D, and for its laters 
tion b — distance of points D to the programmed plane, i.e., plane 


es (Fig. 9). 


Bvare 
) devine 


Tet us introduce brief designations 


2(6) =z, y (3) = ye 2 (0) =z, 
¥,(3) = 9y,, Vy(S)—=vy, 42(9) =P, (2.Tody 


for values of coordinates of the missile and projections of its 


velocity relative to the starting system of coordinates rps ut the 
instant of the cwitening off of the engine ¢ = o. In this :ase it 


MN 
ua 








Fig. 9. 





is considered that the instant t = 0 is referred to the beginning 
of motion of the missile from the start. With fixed values Bs; y, B 
and »,, Vy» ¥, range 2 depends on the width of the place of the _ 
Start and on the location of the plane zy relative t6 the point of 
the compass. The position at the instant of termination of the 
active section of the celestial bodies — moon and sun relative to 
the earth, practicaily does not govern the flight range of the 
missile. A certain small effect on the magnitude of the range of 
the missile proves to oe the state of the atmosphere, especially 
near the points of ivs fali, which subsequently, however, Is riot 
considered. 


Having in mind any specific position of start of the missile 
and its assigned point of fall, we take tha: range 2 is only a funce- 
tion of coordinates x, y, 2 and projections of velocity vis Vis v, 
uf the missile at the end of the power-flight section of flight and 
due.. not depend on the duration of the latter. Thus, we will cone 


Slut that 


Eww LZ, Ys E) Vay Vy, Us) (2.1.2) 


wr 


The construction of functions l(z, y, 3, vie Us v_) in general 
in the form of an explicit function of their variabies is impossible, 


sinee, in tontrast to that given in Chapter i, the eaucetion of motion 


oe et te ae AE ie a ie eat A 





of the missile in the free-flight section in quadratures are not 
integrated. Nevertheless, mathematical machines with great accuracy 
allow comparatively rapidly calculating the magnitude of range of 
the missile 2 according to data of its coordinates x, ¥, 2 and pro- 


jJections of velocity vie vy vl, which corresp.ends to the end of the 


power-flight section in the starting System of coordinates xyz. 

Let us designate by x*, y*, 24 and ve v*, vi those values of 
ccordinates and of projections of velocity of the missile in the 
starting system of coordinates xyz, which corresponds to the end of 
che power-flight section in the calculated motion and by 72% the cal- 
culated range of its flight. In accordance with equation (2.1.2), 
we have 


fea d(s’, Vo 3's Va, Vy, B;)- (2.1. 3)- 


Let us form expression 


Al ae f(z (6), y (¢), 3(¢), Va (t), By (!), ¥e (t)l -- 
* —2(s", ¥, 8» Bay By» Bs). (2.1.4) 


+l can be accepted as the ballistic function mentioned in the 
beginning of this section. Actually, if we turn off the engine at 
any arbitrary instant t, then value of the selected function corre- 
sponding to this instant is directly the error in the flight range 

of the missile. If, however, the engine is turned off at that 

instant t = o at which this function turns into zero, then, naturally, 
there will not be an error in the flight range. 


Because Of the absence of the equation which expresses function 
(2.1.2) in an explicit form in terms of its arguments, the practical 
use of the difference (2.1.4) as a ballistic function is difficult. 
sonsiderably simpler is another ballistic function, which is obtaines 
from the expansion of function 7 into Taylor series about calculated 
values of its arguments, i.e., about the totality of quantities <x’, 
gr, St, vis v4, and vs. According tc expression (2.1.4), retaining 


only smallness of the first order, we have 


i 
i 
i 
} 
a 


iain inn nc 
? putas 


NUMAN 


yen NEON MRNAS NERE 


soem 


SHUM ASRDURRNN 


reo 


ATONE NEAR NAMEN DEN 





Alm (50212 + WO — VIE + (29-2) 4 
+ leat) 0. + 12 — OI ge + tw — 0 £E. (2.1.5) 


Here the partial derivatives ay he ee grr called ballistic 


coeffietents, are taken at values of their arguments, respectively 
equal to x*, y*, 24, vis v* and vie and, consequently, for the 
assigned flight of the rocket are constant quantities. The ballistic 
2verficients can be calculated, specifically, on high-speed computers. 


It is obvious that with a small deviation in the motion of the 
missile from calculated as a ballistic function one can take expres- 


sion 


8) = 21S + WO VISE + (1 FE + 
+172 (t)— vi} zt 17, ()— vy zt TAO oe (2.1.6) 


i 


nctually, with the switching off of the engine at the instant when 
this expression becomes equal to zero, the difference between the 
real and calculated flight range of the missile proves to be a 
magnitude wu: the second order of smallness.’ In this case deviations 
in the parameters of the actual motion of the missile from the cal- 
culated «t the end of the power-flight section of its flight are 


taken az smallness of the first order. . 


§ 2. Equations Connecting Coordinates of the Missile 
in Starting and Nonrotating System of Coordinates 


Yo eelve probicus «ff inertial guidance, the starting system of 


» ordinates represents certain inconveniences. Specifically, the 
determination of coordinates and components of velocities of the 


‘Tf we do not take into account, of course other facts disturb- 
zmg the accuracy of flight of the missile, specifically, the effect 
.f atmospheric sonditions at the end of the power~flight section. 





missile in this system according to readings of newtonmeters or 
integrators of accelerations is complicated by the necessity of 
calculation of translational and coriolis forces of inertia. Let 

us introduce, therefore, the nonrotating system of coordinates Enc 
with the same beginning at the center of the earth as that for the 
Starting system xyz, which rotates together with the earth. Let us 
assume that at the instant of the beginning of motion of the missile 
t = 0 the axes of system of coordinates Ent and xyz respectively 
coincide. 


Let us designate by Urs uy, and Ur projections of angular 
velocity of the earth on corresponding axes of the system of co- 
ordinates &n¢g and by Ui» Ys and U, ~ on axes of the system zryz. 
For an arbitrary instant of time the following equalities are 
valid 


Ug=U, = lU, U,=U,= mv, 
0, =U, = nu, : (2.28.10 
n—- direction cosines of the eartn's axis and, con- 


where 1, ™ 


’ s 


t 


sOpaentliv, vector Y in the system of voordinatcs Eng or, whien ts 
tne game, in the system of coordinates xy. 


After the time t after the start of the missile, tne cvstem cr 
csordinates xyz will be turned counterclockwise at angle 


pa Vt (ee ocd a 


relative to the system €n&, if we observe revetlictn. on the site cr 
the positive direction of the vector of anguiar -elocity sf the 
e2rth U (i.e., on the side of the star Polaris:. 


Together with the current coordinates of the missile x(t), y(tJ 
and z(t) in the starting system of coordinates xyz, let us introduce 
its coordinates §(t), n(t), C(t) in the system ng. Axes of the 
latter, as follows from the aforementioned, do not change their 


rientation relative to directions at fixed stars. 








The table of cosines of angles between the axes of system of 
coordinates €n¢ and xyz! is the following: 


g " " 8 ; 
—£(1—cosg)t+ (1—cos@)mi+ (1 —cosg)nl—! 
. + cos +nsing —maing, 
y (1—cose)im— (1—cosqg)m?+ (4 — cos 9) nm-+- 
—nsing + cos@ + sing, 
5 (1—cos@)in+ (1—cosg)ma— (1—-cos@)n* + 
+ msing —ising + cosg. (2.2.3) 


According to table (2.2.3), one can determine coordinates of 
the missile in the system Eng, 1.e., quantities E(t), n(t), and oft), 
if direction cosines l, m, n, angle @ and coordinates xrft), y(t), 
z‘t) are known. The appropriate equations have the form 


& (2) = (1 — cosg) @ + cosql x(t) + [(1 — cosep)im — 
—ansing] y(t) + [(1 —cosq) In +m sing] z (t), 


a(t) = [(1 Su eaee) mi +nsing] x(t) + [(1 — cosg) m* + 
+cos@l y (t) + [(1 — cosq¢) mn — L sing] z(t), 


Se) = (1 — cosg) al — msing! z(t) + [(1 — cos) mn +. 
_tisingl y () +11 — cosg) n® + cos] z(). (2.2.4) 


The reverse conversion, which expresses coordinates rft). yft) 
and z(t, in terms of E(t), nf(t/ ana C(t), corresponds to the finite 
turn or the system of coordinates Ent relative to the system xyz at 
tie sane angle ¢ but in the opposite direction. Specifically, it 
45 obtained from the previous equations (2.°.4) if we replace in 
them angle « by -% and exchange places cf quintities € and x, n and 


vy, and 21so ¢ and a. As a result we obtain equations 


z(t) = [(1 -- ccc g) BF 4-cneg] & (t) | (Cl — ene g) m4 
+n sing) y (0-1. [1 — cosq) ni — msing) (0), 
y (2) = ((l— cosq) im — n sing] & (t) + [(t — cosq) m*+ 
+cosqin () + ((f — cos) nnt +Isingl S(t), 
s(t) = [(t — cosg) in +m sing] & (0) + [(1 — cosq)mn— 
—Tsing!y (2) + [1 — cosg) n® + cos@l 5 (2), (2.2.5) 


bsee, cp-olfieally, tne fAopendix on page 98. 
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which, of course, can be obtained directly, using table (2.2.3). 


§ 3. Projections of the Velocity of the Missile 
Relative to the Nonrotating 
System of Coordinates 


“ 


Earlier there were introduced designations vift), ve) and 
v(t) for current values of projections on axes z, y and 2 of veloc- 
ity of the missile relative to the starting system of coordinates 
xyz rotating together with the earth. Projections of velocity of 
the missile in nonrotating system of coordinates &n¢ on the axis of 
this system are designated, respectively, by up(t), u(t) and u(t). 
It is obvious that 


ue SO, (ye BO, n= Sh, (2.3.2) 
in the same way as 


v(t) = $8, 9, (t) = S910. , v,(t) = a0 ° (2. 3.2) 


Quantities v(t), v(t) ana v(t) can be expressed by up lt), ult), 
uy f(t) and inversely. According to equations of kinematics relative 


to the motion of the points, we have 


By (t) = Dg (2) + Uys (t) — Uy (0). 
8, (t) oy (t) + U2 0) — O42 (8), 
w, (6) = v, (t) + Uy (8) — Ux (8), (2.3.3) 


where u(t), u(t), u,(t) — projections on the axes x, y and 2 of 
the velocity of the missile in the rotating system of coordinates 


eng. 


In equations (2.3.3), in accordance with equalities (Zeted is 
let us replace Un uy, and Us» and coordinates z(t), y(t) and z(t) 
will be expressed by means of relations (2.2.5) in terms of co- 
ordinates &(t), nf(t) and c(t). If-further projections un(tss u(t) 
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and u,{t? are expressed, using table (2.2.3), in terms of projections 
ae u_.?t; and u(t), then we will arrive at the following desired 


equations 


pelt) = [(1 — cos@) + cos@] ug (t) -+ 
+ (1 — cose) im +n sing) uy (0) + : 
+((1 — cosg) in — m sing] uz (t) — Uy (I(t — cos q)ni 4 : 
+maingl & (t) + ((1 — cose) mn — Ising! y (0) + 
+ [(1 — cos q) a*+cos ql € (0)} 4+ Us (11 — cosq) im — 
— wsings & (t) + ((L — cosg) m*4.cos@l n (t) + 
+ ((1 — cos@) ma + 1 sing! § (6), 


oft) = {(1 — cos@) im — n sing) ug (t) + 

+ I(1 — cosg) m* + cong] uy (t) + : 

+ [(1 — cose) ma + 2 sing! uc (t)— ° 

— OU, (Id — cosg) & + cose] & (t) + 
+ ((1 — cosg) mi +n sing) n (t) + - ; 
+ (1 — cong) al — msingl § ()} + 
+ Uc{{(t — cosg) in + m sing! & (t) + ((1 —cosg) mn— 
-- Using] y (t) + (C1 — cong) n* + cos@! f (4), 


p(t) = (1 -- cosg) in +m sing] ug (t) + 
+ If — cong) ma — Laing! uy () + [1 — cosg) n® + 
+ coog] ue (t) —- Ue (I(1 — cosq) im ~ n sing] & (t)+ 
+ [{4 — cosg) m® + cose! m (t) + 
+ ((f — cong) mn + 2 sing! ¢ ()} + & 
+9, (11 —cosg) P + cos@) £ (t).+ ‘ 
+ [(1 — cosg) mi +A sing) y(t) + 
+1(1 — cosg) ni — m sing) & (¢)). (2.3.4) 


The last equaticns are somewhat simplified after the repiacement in 
then of projecticns of angular velocity of th earth Urs uy and Uy 


respectively by products lU, mu and nU in acccrdance with the same 
equaiities (2.2.1). As a result we obtain 


og(t) = (1 — cong) F + cos) uy (t) + 
+ [(1 — cosg@) mi +n sing] we (t) + 
+ ((1 — cong) nl — m sing] u, (t)— 
— U {coag Img (t) — an ()1 — 
— basing UE (0) + mn (6) +06 (t)) + & (1) sing), 











DB, (t)= I(f — eosq) im — nsing! ue (0) + 
+ [(1 — cosq) m* + cosg] uy (t) + 
+ (1 — cosg) mn + lL sing] ug ()— 
— U {cosg Ink (t) — 10 ()) — m sing lE () + mn (2) + 
+ n§ (t)] +7 ( sing}, 


o,t)= [(1 — cosg) In +m sing] u(t) + 
+ [(1 — cosq) mn — I sing] uy (t) + 
+ [(1 — cos) n® + cos) ug (t) — 
— U (cos [im (t) — mE (0) — 
— nsing (28 (¢) + mm () +8 (¢)1 + € (¢) sing). (2.3.5) 


It is obvious that by similar means it is possible to arrive at 
equations which express quantities Ups uy and ue in verms ofr vi 


v «ind vo. we nave 
y 2 


te (t)= [(1 — cosq) # + cose] vy (t) + 
+ [1 — cos@) im — a sing] v, (t) + 
+ (1 — cosq) In + ming] », (t) + 
+ U {cos@ [mz (t) — ry (t)} + 2 sing [tz (t) + - 
+ my (t) +172 (or— z (t) sin 9}, 
uy (t)= [(1 — cos) mi + nsing] vz (t) + 
+ [(1 — cosq) m* + cos@] vy (t) + 
+ [4 —cosg) mn—tsing) v, (t) +7 (cos pinz (t) ~ 
— Uz ()] +m sing [Iz (t) + my (t) +22 (I—y sing), 


uc (t)= [(1 — cosg) nl —m sing] v, (t) +4- 
+ (1 — cosg) nm }'Esingl v, (0) + 
+ [(1 — cose) n® + cosql », (t) + 
+ U (cosq [ly (t) ~— mz (t)] + 
+n sing [iz (t) + my (t) +z ()] — z (4) sing}. cds 


5 4. Ekrror in Flight kange of the Wissile #3 4 
runetion of Changes in Furameters of the Ena 
of the FPowered-Plight Seation in the 
Nonrotating System of Coordinates. 
Initial Ballistic Function 
We will designate coordinates and projections of velocities of 
the missile in the nonrotating system of coordinates €nt at the 
instant of the end of the powered-flight section of its flight 
= ', respectively, by letters €, nh, & ama w,, u, gn aes Die 


—_— 


are similar to designations (2.1.1) introduced earlier for quantities 
2, y, 2 and vee Pye v,» which refer to the starting system cyz. Thus, 


E(o)=% nan, (9) =f, 
wy (0) me te, Hy (0) me try, tg (8) am te, (2.4.1) 


Setting in equations (2.2.5) and (2.3.4) ¢ = Uo, we obtain expres- 
sions of quantities x, y, 2 in terms of —, n, © and o and also 
quantities v., v, v, in terms of Ups Uns Urs Es, nn, G ando. This 
allows considering magnitude of range 2 as a function of values of 
coordinates of the missile &(t), n(t) and f(t) at the instant t = o. 
projections of its velocity up(t), u(t) and u(t) in the nonrotational 
system of coordinates EngZ (at this instant of time) and, finally, the 
duration itse!* of the powered-flight section o. In accordance with 
formula (2.1.2) we now have 


len l(z,y, £, Dg, Dy. 0,) = 1 (E, The Or Ue Uns Ue’ 9). (2.4.2) 


it is easy to explain why in this case the range 2 clearly depends on 

Variable 9. Actually, the position of the rocket relative to the 

earth at the same values of its coordinates in the nonrotating system 

Eng sucst2ntially depenas on the position of the latter with respect 

to starting system of coordinates zyz connected to the earth, i.e., 

on angie ¢ = lo. The came refers to magnitudes of projections of i 
roof velocity of the missile in the etarting system at the 


t‘netant of termination of the powered-fiighi section. 


Calculated values of parameters of the end of the powered- 
Tlight section in the starting system of coordinates xr*, w*, 2%, 
v*, v*, and v* correspond to calculation values of parameters g*, 
n*, ce, urs “a, and us in the nonrotating system Ent. They are 
connected with each other by equations (2.2.5) and (2.3.4) or (2.2.4) 
and (2.3.6 , in which it follows to assume ¢ = Uc*, considering o* 
as the Jesignation of the calculated duration of the powered-flight 
secticr of flight of tne missile. It is ubvious that in uccordance 


. t. 8, aba 4 : = 
situs ruUPmule £ oie Sy 


has 
Wy) 


2, A Aka. ote 8 


a8,» eng 





Bm l(sy', 2, 05, 0p Bs) = LE WO, ues ae 4es OD, (2.4.3) 


where 2*, as previously, is the rated value of the range. 


The difference 4Z between the magnitude of the actual range of 
the missile 7 and its calculated value 72* can be represented, similar 
to equation (2.1.5), in the form of the expansion 


A= E—EV t+ (M— MI et (C—O) H+ (we — 0d) Se + 
Ht HF + (te — 0) FH + O— ODF, (2.4.4) 


* 


where terms of the second and higher order relative to differences 
are dropped 


t— e, a—", t—e, ug — Be, a a, | ad Re, o—’, 


In contrast to the expansion (2.1.5), entering into equation 
(2.4.4) is a term of the first order, which contains by a factor 
the difference o - of between the actual time of the powered-flight 
section of flight of the missile o and its calculated value o*. The 
same difference is contained by a number of terms of expansion 
(2.4.4) of a higher order. The reason for this fact was already 
explained in the beginning of this section. 


. & #& & MB « a a, 
Derivatives Xi Fa! Hi Im’ Sans ing and 3 are themselves functions 


of variables —, n, f,; Ups Uns My and o. In the expansion of (2.4.4) 
they should be taken at calculated values of enumerated variabies 

ER Ht Bh, ugs uss ues o*, and, consequently, for the specific 
assigned flight of the missile are constant quantities — ballisti 
coefficient, which refer to the nonrotating system of coordinates. 
Let us show that they are all expressed in terms of ballistic 
coefficients of the starting system * =. ie i and i 
Actually, according to rules of differentiation of complex functions, 
specifically, we have: 


a3 








a at as , ata a 2a a are : 

RBA HHH tet e 

at as Oe al ov at Oe, 

Bay ~ 50, Sug + Be Bu + 30, Bue Senn) 


Analogous expressions take place for derivatives ee Bat bag’ 


and for the partial derivative = equation 


at 4M bs 8! by a as, at Or, , at ae, , Ot bv, 


1s valid. Entering into enumerated expressions, the partial deriva- 
tives cr variables xr, y, 2, V2 ay v, with respect to variables €, 
TN, &s Ups Uns vy and o are formed by means of the application of 
appropriate operations of differential equations of the type (2.2.5) 


and (2.3.5), in which it follows preliminarily to set t = o and 
¢ = Uo. 

Thus, for instance, 
ov, 


Hw Bes = C8 + (4 ~ cos Us) 2, 
a, 


™ 5, = (1 — cos Vo) im +nsinUs, 


* ae ‘ay 


=O (t—P)sin Vo, 


= we U ([—(1—P)EsinUs + (im sinUs +-neos Us) y + 
+ (in sin Us — moos U0)E}j, 

z wa T[—(1~ 2) ugsin Us+(lm sins -- ncos Ua) 24 + 
+ (Insin 0s — moos Ua) ug] + U* (mG — an) sin Us + 


+E + my + nf) ers Us — Ecos is} C204. 


if now in equations (2.4.5), (2.4.6), (2.4.7) and ones similar 
to them, we assume the values of ali variables equal to the cal- 


cuiated values, then thus the desired coefficients in expression 


(2.4.4) will be. completely determined. 


8 4 





We will assume that as a result cot’ the small deviation in the 


real motion of the missile from the calculated in the expansion of 
(2.4.4), the terms caused by the so-called lateral deviation in the 


beg eh eaten 8 = at sy 2 * 
missile @— 1x and (ug — ue) 5 are smail in comparison with other 


Sorms. \S a result the difference between the actual und calcu!uted 
Vilght range of the missile Ai can be represented by the apprcximate 


syuation 


Al = &—*) 5 +(n—nIS + (me — ay) 2 + 


a ; 
+ (ty— 04) 5. +(0— > 2.4.8; 

oe at oat at f 
where su, and gp» aS was already shown above, are numbers 


a dug’ du, 
which are completely defined by the selected calcuiated motion cr 


eee 


to the foilowing owviiisii. 


¢ 


& 
‘wrsstiun, whieh we call the inttial ballistic function: 


a(t) = nO >iia + In) — WISE + tne () — vel Set 


+a — ui +00). ai 


Sweytee fed fo de 
ee DS? Seu Se lg ", rw 


» -(ét) and un(t), us nefore, are movir. evecorainutes 
re 


md projections of velocity or the missile in the nonrotating syvctem 


DY osGurdinautes Eng. 


If the engine of the missiie is turned off ut any cn.0tunt ° = ty, 
‘tisthy. “G, the ater er equrtirne (peli) and (eat) the sae ota as 
mGddistie functton (2.4.9) corresponding to this Instant correct ot. 
osmabainess of tue second order determines the error in tre filent 
range of tne missiie. To avold this error, one snould turn off tne 
engine of the missile at the instant when tne tallistis runcticn™ 


‘oo *urnms into sero. Thus, tne necessary instant of the switening 
3 od 


is determined by tne root of equation 


weet he tietine is mr 








RO-ENZ + nO — 1K + ae — 0 ZS + 


+e) —0 Zo) =O. (2.4.10) 


Equation (2.4.10) will be called the initial ballistic equation. 


The basic content of subsequent sections of this chapter anid 


‘also Chapters III and IV consists in the construction of a series 


of 


as 


‘the flight 


other ballistic functions, with the help of which it is possible 
simply as possible to solve the problem about the inertial control 
range of the ballistic missiles. They all are obtainec 
means of certain conversions from function (2.4.9) and different 


presentations of coordinates and projections of velocity of the 


ef 


Tek 


Ce 


moving |. 
Terential 
readings 
snould 


fF and mn of the nonrotating system of courdir :tes Eng. 


missiie in the nonrotating system €n¢ in terms of current readings 


of 


integrators of accelerations. 


§ 5. Differential Equations Which Determine Current 


Cocrdinates of Motion of the Missile 


The most natural, but, as it appears, not the simplest, methcd 
eonstruction or the ieft side of the initial bailistic equation 


(2.4.10) aboard the missitie is the preliminary obtaining of most 


é 


‘ourdinates &(t) and n(t) by means of solution of the dif- 
equations, which connect these coordinates with current 
two newtonmeters. The axes of sensitivity of the latter 
fixed directions, respectively, parallel to axes 
Readings of 


ofr 
maintain the 


the mentioned newtonmeters are projections ¢.¢t) and a(t) of the 
S : 


apparent acceleration of the missile on the axis € and n. According 
to equation (1.4.1) of Chapter I, in this case we have 
= 7 _ tj) = t)— . ‘ 
facile Gz (t) = wy (t)— fe, a,(t) tent) fe (2.5.1) 
Wrins ote . | 
we (t) = A, w() = SO (2.5.2) 
ire axes € and n of the real acceleration or the 


t+ the flxed system cf coordinates Sn¢é, and 


Was 
nm 








fe lB), 00), B05 41. 
fam falB(, 100), (0 1 


2 


(2.5.3) 
; 


\ 


are, respectively, projections of acceleration of the force of the 
earth's gravity, which depend in evident form not only on coordinate 
of the missile, but also on time ¢. For an explanation of the last 
fact let us note that even at the point with constant coordinates 
relative to the nonrotating system Ent the acceleration of the force 
of gravity is changed with the course of time. Actually, with 
respect to the system of coordinates Eng the position of the earth 
is continuously changed. At the same time the field of gravity of 
the earth does not have radial symmetry, specifically, as a result 
of the deviation in its form from a sphere. If the last fact is 
disregarded and the earth is considered a sphere with radial and 
symmetric distribution of density, then projections of acceleration 
of the force of the earth's gravity will be represented by equations 


ha—pRiD, pe — pad, (2.5.4) 


no longer dependent in evident form on time, Into them fo is the 
value of acceleration of the force of gravity on the earth's surface 
examined as a sphere of radius R, and 


p= VEO MOP + emi (2.5.5) 


is the distance between the missile and center of the earth. Becauss 
of the smallness of the lateral deviation in the rocket %/t in com- 
parison with the radius of the earth R, this coordinate in equation 
(2.5.5) can be cropped. Correspondingly, it is possible to drop 
coordinate &(t) in equations (2.5.3), considering the effect of this 
coordinate on the magnitude of acceleration of the force of gravity 
unimportant. Then relations (2.5.1) can be considered as a combina- 
tion of two differential equations 


FBP am on (6) + fe 1B WED: AI. 


' | (2.5.6) 
FAO wo ag(t) + falb CE) ACO: 
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for the search of destred functions &(t) and nf(t/ on the given 
current readings a(t) and aft) of two newtonmeters. 


The direct colution of the system of nonlinear differential 
equations (2.5.6) aboard the missile for the purpose of seeking the 
current values of coordinates €(t) and n(t) it is possible by means 
of the use of a rather high-speed computer. In the subsequent two 
chapters the approximation solution to these equations, which in 
principle can be used in the system of control of the flight range 
of the ballistic missile in the presence of onboard of its newton- 
meters or integrators of accelerations and simplest computers is 


used. 


§ 6. Auxiliary Relation Connecting Magnitudes of 


Ballistic Coefficients in the Nonrotating 
System of Coordinates 


At different conversions of the initial ballistic equation 
(2.4.10) one relation proves to be useful, and it connects magnitudes 


of the ballistic coefficients it He im i and 2 in the nonrotating 


system cf coordinates Eng with projections of velocity w#, un and 

with projections of acceleration of the force of gravity f, (ee, n*; 
o*), £,f8" n*; of), which refer to the instant t = o* of the termina- 
tion of the powered-flight section of the calculated motion of the 
missile. 


Let us introduce functions &*(t) and n*(t), which are current 
coordinates of the missile in its calculated notion on the powered- 
flight section. They satisfy the totality of differential equations 
(2.5.6), if in them the current values of projections a, (t) and a(t) 
of the apparent acceleration are replaced by their calculated values 


a 2 a 
agft) and antt). Thus, 


a(t) oe : ae 
Sie (t) + A 1B Ct), 1° (4); 41, 
" . (2.6.1) 


Se = a3 (t) + FIED, (0; A 
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Correspondingly, let us introduce functions &**(t) and n**(t), which 
are coordinates of the missile in the free-flight section of its 
calculated motion. Let us assume that in the beginning of motion 

of the missile in the free-flight section the resistance of the 
atmosphere is unimportant, and, consequently, the apparent accelera- 
tion is absent. In this case one should consider that functions 
E**(t) and n**(t) satisfy the totality of differential equations, 
which is obtained from (2.6.1) if in equations of the latter we drop 
quantities ag(t) and an(t) and replace variables €*(t) and n*(t) 
respectively by &€**(t) and n**(t). We obtain 


Seem felB" (0), 7 (0 EI. 


., (2.6.2) 
SHO mm fa18" (0, 9 OE 


Let us designate the current values of projections on the axes € and 
n of the velocity of the missile in the powered-flight section with 

its calculated motion by ug(t) and un(t)s and by up*(t) and unt(t) - 
the corresponding quantities which refer to the free~flight section. 
We have the evident equalities 


ne SO, k= SO (2.6.3) 


and, similar to them 
ng (t) = 8, us (t) = SO, (2.6.4) 


The powered-flight section of the calculated motion of the missile at 
the instant t = o# will pass over into the free-flight section, whence 
it follows that 


B°(0*) = 8°" (0°) = E°, 1° (3°) me H* (6") me 1°, (2.6.5) 


and also 


ue (3") = ug (0°) = ur, ty (3°) = uy (3°) ae ay. (2.6.6) 





According to equations (2.6.1) and equations (2.6.5) for the end of 
the caleulatei powered-flignt section,'! we obtain 


FES) at (0") + fe(B ns"), 
o.. a : (2.6.7) 
(as ) ro a,(3') + hk (, "3 a). 


Similar relations can be ootained on the basis of equations 
(2.6.2) and the same equations (2.6.5) for the beginning of the 
calculated free-flight section. They have the following form: 


SoA) ae ee, nv; 3’), ae. W(B, wi 3°). (2. 6. 8) 


Let us examine now on the calculated free-flight section of 
flight of the missile the certain instant ¢t = tas which follows 
after the instant of the switching off of the engine t s o*. Correct 
te srallness of the second order we have expansions in Taylor series 
fer coordinates of the missile fy and Ny which refer to this instant, 
namely: 


Ber y= to) + (4-0) SE, 


mm A(t) en + (9) SS), eae 
or, taking account equalities (2.6.4), (2.6.5) and (2.6.6), 
Ra + (4s— oa, mee +(— ou. (2.6.10) 
Similarly for projections ue and ui of the calculated velocity of 
tre missile av the mentioned instant t = t, we have 


1 


‘in designations of the type d&(c)/dt it follows, of course, to 


consider, that from the beginning an operation of differentiation of 
the aperopriate function is produced, and then the vaiue shown in 
parenthesis ic given to the argument. 


e407 pm pane Ya GN oe TMG AD EMEP =e 8 





ome PU CUTTORE URAeeslmincat TE 








t= nTe= ao +,—o LE, 
wt mals) marin) +o). 


(2.6.1) 


On the basis of equations (2.6.4), it is possible here to produce 
the following replacement: 


SO A OL A, (2.6.12) 


Taking into account, furthermore, equality (2.6.6), and also rela- 
tions (2.6.8), we obtain 


Bt — B+ (la — 9) feb, WO), 


(2.6.13) 
Bm Bt (ta — oD FE. WS). 


Let us imagine such a motion of the missile in che powered- 
flight section different from the calculated, as a result of which 
at the instant ¢ «= ty it proves to be at the point with coordinates 
gs and No possessing a velocity the projections of which on axes 
€ and n are respectively equal to quantities Ke and es if at this 
instant t s ty the engine of the missile is turned off, then the 
following free-flight section of its motion, because of the unique- 
ness of soluticn to the problem of dynamics completely coincides 
with the calculated. Actually, both motions are such that at the 
same instant t = ty in them positions and velocities of the missile 
relative to the system of coordinates Ent coincide. Purthermore, 
when t > ty they are subordinated to the same totality of differential 
equations, which refers to the free-flight section of its motion. 

p Thus, at such a powered-flight section the rocket will not nave an 
error in the range of its flight. Consequently, parameters of the 
end of the powered-flight section of its motion €)> Ny Me and ur 


should satisfy the initial ballistic equation (2.4.10), i.e., 


RN +a— oF + elu H+ 
+E + (h—o ZS m0. (2.6.24) 











Replacing here quantities By and ny by their expressions according 


to equations (2.6.10), and ue and ue in accordance with equations 


(2.6.13), we obtain after, reduction to the common factor t.- o* 
equality 
BTA TAG TONE + 
+AU. Wi OLE +E = 0, (2.6.15) 


which connects the ballistic coefficients x x. Fes bes and Bad in the 


nonrotating system of coordinates €n% with projections of accelera- 
tion of the force of gravity (8. 95 0°) fs (8. na and with projec- 
tions of velocity of the missile ug and ur at the calcuiated instant 
cf the switching off of the engine. 


Let us exclude from equality (2.6.15) projections of accelera- 
tion of the force of gravity by means of relations (2.6.7) and use 
eguations (2.6.3). As a resuit let us obtain the new equality 


Be) Mg SI) Bg (PE) car] 
+[P9 Ed —an(on] FZ + =O. (2.6.16) 


which will be used in the following section during the conversion of 
the expression of the initial ballistic equation (2.4.10) te the form 
convenient for applications. 


§ 7. Isochronal Variations of Coordinates and 


Projections of the Velocity of the Missile. 
Basic Ballistic Equation 


With the juxtaposition of the real and calculated motions of 
the missiles, let us call the isochronal variations of its co- 


ordinates the differences 


BE (t) =m E(t) — F(t), 8n(¢) = H(t) —0°(¢). (2.7.1) 
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Similarly, let us introduce isochronal variations of projections of 
velocity of the missile 


Sug (t) = ae (f)— E(t), dun (t) = By (t) U5 (2) (2.7.2) 
and projections of the apparent acceleration 


Bag = Gg (t)— ah (t), dg = a(t) — 05,(1). (25/103) 


We will consider the enumerated isochronal variations small quantities, 
the squares and products of which can be neglected. Thereby, it is 
proposed that the problem of thrust of the engine both in direction 
‘and.in magnitude is produced to a sufficient degree accurately. 


Functions €*(t), n*(t), ug(t), un(t), agit) and an(t) refer to 
the powered-flight section of the calculated motion of the missile 
and, consequently, are determined at values of the argument ¢ not 
exceeding the duration of this section o*. At the same time in 
equations (2.7.1), (2.7.2) and (2.7.3) functions €(t), n(t), up (ti, 
u nft)s a(t) and ay, (t) refer to the powered-flight section of the 
real motion of the tea iies the duration o of which can be both less 
and more than the calculated value o#. Therefore, for the complete 
certainty of variations s&(t), é6n(t), Sup (tl, buy, (t), Sap (t) and 
Sa, Che in the whole time interval of the enaanea. flight section of 
the real motion of the missile, one should agree upon what is under- 
stood by functions €*/t), n*(t), uslt), us(t)s az(t), an(t) with 
argument t somewhat exceeding quantity o* (within limits of the 
allowed variance of duration of the powercd-flight section for t 
specific type of missile). If one assumes that the thruct ef the 
englne with termination of the powered-fiight secticn is lcwerea 
eridully or by stages, then for these functions when = > OR POLS 
natural to take, respectively, functions &**(t), n**(t), ug*(t) and 
unr(t), which refer to the calculated free-flight section and assume 
eineetous ag(t) and an(t) equal to zero. In this case it follows to 
expect that the tsochronal variations (2.7.1), (2.7.2) and (2.7.3) 
will be small at the short time interval directly following the 
calculated Instant of the switching off of the engine. 
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For missiles with the sudden switching off of the engine, on 
the contrary, it is more preferable to consider functions &*/(t), 
n*(t), ut(t) and un lt) when ¢ > o* equal to their calculated values 
calculated on the assumption that the. engine of the missile at the 
instant t = o* was not switched off. Here the isochronal variations 
(2.7.1), (2.7.2) and (2.7.3) again can be considered small. 


Let us use now equations (2.7.1) and (2.7.2) for isochronal 
variations of coordinates and projections of the velocity in order 
to transferm the initial ballistic equation (2.4.10) by which, as 
was already mentioned, the necessary time of cessation of operation 
of the engine of the missile for its hitting of the assigned target 


is determined. 


Let us present in the beginning equation (2.4.10) in the 


following manner: 


ROB OSE t IO) — VO +0) — HO 
+i = — [OES + 
Hin a1 + OME + 
+192 Re +E F. (2.7.4) 


side of the last equation let us substitute expressions 
variations cf ccerdinates and projections of the 
ket, according to eguations (2.6.5), (2.7.1) and 


Into the left 


velocity oi C dr 
Loe Fee? In its right side lct us cxpand each of the functions §*(t), 
merits, weitl, ance in Taylor sertcs near the value t = s*, retain- 
$yg in the cxpansion only terms cf the first order of sirallness. we 
obtain 


83S + OnE + bu; i + au, () 2 = 


~ aun Es 
reece ee (2.7.5) 


+ aa ou, 


The right side of the last equation can be considerably 


simplified if we take into account equality (2.6.16) of the 
previous section. 


As a result we arrive at the following equa- 
tion: 


BRE) Se + Onl Ze + Ove) FE + Oa Gem 
= -0 [QO +00 FE]; (2.7.6) 


which, sudsequently, we will call “he baste ballistic equation, 


For the actual use of equation (2.7.6) in the system of inertial 
control of the flight range of ballistic missiles its left side should 
be constructed aboard the missile at current readings of newtonmeters 
in the form of a certain electrical or mechanical magnitude. 














CHAPTER Iil 


APPROXIMATE PRESENTATIONS OF THE BA'ULISTIC EQUATION 
DETERMINING THE INSTANT OF SWITCHING 
OFF OF THE ENGINE 


§ 1. Differential Equations for Isochronal 
Variations of Coordinates of the 


Missile in the Powered-Flight 
Section of tts Motion 


In the previous chapter there was obtained tne so-called basic 
bialistic equation (2.7.6), which snould satisfy parameters of the 
end of the powered~flight section in order that the rocket would 
not have an error in the range of its flight. For the construction 
of the left side of this equation, it is necessary to know aboard 
the missile the isochronal varia:ions of coordinates 6€(t) and én(t),, 
ana alse projections of the veiocity Su, (t) and bu, (t). The latter 
are, of course, time derivatives of variations of coordinates, i.e., 


=. 


t=O, ta) a BY. Geld) y 


The coordinates themselves E(t) and n(t) «f the real motion of 
the missile in the powered~flight section satisfy the totality of 
difterential equatilun. (7.5.6), and coordinates of the calculated 
Motion — totality (2.6.1). Let us form the differences, respectively 
of the left and right sides of these equations and equate them to 
each other. Taking into account equalities (2.7.1) and (2.7.3), we 


= 


obtain, 
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FEO aut) + felB UO + 88s ME + On; — 
ALO WO: A 


POT) me Boat) + fa 1B (6) + 8B (), (4) + ON; HI — 
— Fal), 8 (0; A . esis 


In the right sides of the last equalities we can use the 
expansion in the Taylor series of the function according to increases 
in two of its arguments. Retaining only terms of expansion of the 
first order relative to variations 6&(t) and é6n(t) and carrying out 
simplifications, we will arrive at the following totality of two 
linear differential equations: 


PEO boy (t) + ME ag ce) + FE ancn, (3.1.3) 


oo ox 8a, (t) + mf dE (t)' + Fn. 


The desired functions of these equaticns are isochronal varia- 
tions of coordinates of the rocket 6&(t) and 6n(t). Initial condi- 
tions for them can be obtained from the fact that at the instant of 
the launch, 1.e., when t » 0, coordinates and velocities of the 
missile are the same both in real and in calculated motions. Con- 
sequently, according to equations (2.7.1) and (2.7.2), we have 





88 (0) = 0, $n (0) = 0 (3.1.4) 
and further 
bug (0) = EO 20, 340) = 9) 0. (3.1.5) 


The right sides of equations (3.1.3) contain variables — 
isochronal variations of projections of the apparent acceleration 
of the rocket Sa, (t) and ba (t). In accordance with equations 
(2.7.3), in principle they can be obtained aboard the missile by 
means of a continuous formation of the difference between the real 
reading of the appropriate newtonmcter and its current calculated 


vaiue. 


KT 





Coefficlents of equations (3.1.3)? 


e Bah at (3.1.6) 


should be considered as time functions, known earlier for each 
specific case of flight of the missile. In order to obtain these 
coefficients, in accordance with rules of expansion in Taylor series, 
in partial derivatives of functions f# (89; 4) and f.(& 034) according to 
variables — and n, it is necessary to produce replacement of the 
latter by functions &*(t) and n*(t). Functions &*(t) and n*(t) are 
earlier known according to calculated motions of the missile in the 
powered-flight section. Thus, if aboard the missile continuous 
(without substantial lag) integration of equations (3.1.3) is carried 
out and thus magnitudes of functions 6&(t), én(t), Sup (t) and bu, (t) 
become known, then the left side of the basic ballistic equation 
(2.7.6) can be constructed by means of multiplying and adding devices. 
Similcrly, with the use of a clock mechanism the right side of this 
equation is constructed. 


§ 2. Approximate Expressions for Isochronal 
Variations of Coordinates ana Projections 
of Velocity of the Missile. 
Sinpititeation of the 


allistic Equation 


The simplest method of the use of integrators of accelerations 
in the system of control of the flight range of ballistic missiles 
is based on the simplification of differertial equations (3.1.3). 
In them terms containing as factors the isochronal variations them- 
selves of coordinates 6&(t) and 6n(t) are dropped. This is equivalent 
to the assumption about the fact that the effect on the missile of 
forces of gravity with its actual and calculated motions can be 
“sonsidered practically equal. Such an assumption is admissible only 


‘The equality in (3.1.6) follows from the fact that projections 
cr acceleration of the force of gravity are partial derivatives of 
tne potential of gravity according to appropriate coordinates. 


ke 


with sufficiently small deviations in the real motion of the missile 
from the calculated, i.e., at increased requirements for the accuracy 
of metion cf the missile according to the assigned program. Equations 
(3.1.3) are replaced in this case by the approximate equalities 


SHO mn du(, — SO = 80,00), (3.2.1) Here 





path 
whence allowing for relations (3.1.1), (1.5.4) and (3.1.5), there are 
equations for variations of projections of the velocity of the missile 
where 
, e 
Sag (¢) ao m \beccedat = OV e(t), 
(3.2.2) 
a 
ee (t) = SO = Vda, dt = aV a(n agecn 
Quantities SV, (t) and éV Ct), which enter into these equations, 
can te c&llieu isochronal variations of tne upparent velocity cf the 
. missile. According to equations (2.7.3), they are differences thet 
i Belt) = Ve(t)— Me), SV a(t) = Malt) — Vac, (3.2.3) 
tivel 
‘ in which 2 ; equat 
Vata faint, Vat Sand Bok of th 
‘ , 
are projections of the apparent velocity of the missile during its 
actual motion and 
} ‘ . s 
Fe) = Jaca VAQ= Jescna (3.2.5) 
| and a 
and ( 


are calculated vaiues of the same magnitudes. 


Integrating in turn expressions (3.2.2) with respect to time, 
we obtain the fotlowing equations for isochronal variations of 


ccorduinates of the missile: 








te t 
88 (t) = \} Bag (1) det Javea = BS5¢(1), 
oe 1 (3.2.6) 
an(t) = \ be, (t) di? = Savana = 85,(!). 


Here 6S, (t) and $5, (t) are variations of projections of the apparent 
path, i.e., differences 


BH = HHO— SKN, SAW = Sa(— SA), (3.2.7) 
where , 


: : ie 
S(t) = \vecnar, S()= J vaca (3.2.8) 


are profsections of the apparent path of the missile and 


So) = friar, Saft) = \ vee (3.2.9) 
their calculated values. 


Let us substitute now expressions (3.2.6) and (3.2.2), respec- 
tively, for 6&(t), én(t), Sup (t) and Su, (t) into the basic ballistic 
equation (2.7.6), which determines the instant of the switching off 
of the engine. We have 


OS  + OST + ON + OV ONE. 
= [QOS +075] (3.2.10) 


Replacing here variations of projections of the apparent velocity 
and apparent path by their expressions, according to equations (3.2.3) 
and (3.2.7), we obtain 


1S) — SIZ + 15 — SO Fe 
° a . 
+e —VEOLR, + (Ka) — VA) ie - 
= —¢ oy [orgy + 2,672]: (3.2.11) 
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The last equation can be considerably simplified, if we note 
that correct to smallness of the first order inclusively the follow- 
ing expansions in Taylor series take place: 


Vell) = Vile’) + (t—3') ac (3°), 


oye e, eo @ e. ( 3.2. 12) 
Va) = Vio) + (¢—0') a, (9°). 


By means of these equalities equation (3.2.11) is reduced to the form 


1Se () — Se)! a+ [Sa(t)— (9) x +Ve(t) a + 


+ V0 = VEC) + Vale) S. (352033) 


where terms dependent on time in evident form are already absent. 


Thus, the basic ballistic equation (2.7.6) can be substituted 
by the approximate equation (3.2.10) or the equetion (3.2.13) 
equivalent to it. For the construction of the left side of equation 
(3.2.13), the presence is necessary aboard the missile of two 
integrators of accelerations, axes of which during the whole powered- 
flight section of flight should retain directions, parallel respec- 
tively to axes € and n of the nonrotating system of coordinates &nc¢. 
Furthermore, there must be a computer, wnich includes in its com- 
position two additional elements for the time integration of current 
readings of the very integrators of accelerations and a special 
element for the reproduction of the calculated current values of 


projections of the apparent path Set) and Sift. 


§ 3. Construction of the Ballistic Equation by Means 
of the Use of Readings of Two Integrators 


of Accelerations with Special 
Orientation of Their Axes 


of Sensitivity 
There can be cor.siderable interest in the possibility of the 


reduction in the number of elements of the computer system of 


inertial determination of the flight range of the ballistic missile 
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by ineans of the proper selection of orientation of axes of sensitivity 
of integrators of accelerations. It turns out that there can be 
produced a total of one additional integration of current readings 

of integrators of accelerations, and not -two, as in the method of 

the construction of the ballistic equation given in the previous 
paragraph. For this purpose, as will be shown below, axes of 
sensitivity of integrators of accelerations should be parallel te 
certain directions fixed in each specific case of flight of the 
missile (so-called A- and p-directions). 


Another method of simplification of the system of inertial 
control of the range with the use no longer cf an integrator of 
accelerations, and the meter of the apparent accelerations with the 
changing orientation of the axis of sensitivity, is stated in the 
following paragraph. Subsequent modificaticn of this method is 
Ziven in &§ 5 of this chapter. 


Tuus, let us position the axis of sensitivity of one of the 
integrators of accelerations in the plane &n of the nonrotating 
System of coordinates €nt at a certain constant angle A to the axis 
€. The direction itself of the axis of sensitivity of this integra- 
tor will be called the A-direction (Fig. 10). 


On the basis of equation (1.6.5), after the replacement in it 


of letters x, y and v, respectively, by €, n and ’ and also angie a 


v 


by angle 4, we obtain equation 
Vy (t) = Vg (t)cosd + V,(¢)sind, (3.3.1) 


whith expresses the projection of the apparent velocity of the 
missile on the mentioned A-direction by its projections on axes € 


and mn. It 1s obvious that in accordance with formula (1.4.7) 
’ 
Vi = Facade, (3.3.2) 
where a, {t) — component of apparent acceleraticn along tne same A- 
wdireetion, 





Fig. 10. 


Thus, quantity Vy (t) is the current reading of the integratcr or 
accelerations with the axis of sensitivity, oriented in directicn ). 


Let us select angle ’ so that there would take place equalities 


Fim Moos, Roo Mind. (35353) 


ere, as it is easy to see (see Fig. 10), quantity 


Ka Fe) +(s.) (3.304) 


is for each specific case of flight of the missile, an earlier 
Known quantity. 


Using equations (3.3.1) and (3.3.3), let us transform the sum 
of tne last two corresponding terms of the left side of equation 
(3.2.13). We have 


"Far + Vee wo M [Ve (t)eosd + Va (tain A] = 
o MV, (0 cu oe5) 





Tnus, the reduced sum correct to the constant factor M is 
determined by the reading of the integrator of accelerations, the 








axis of sensitivity of which is located along the A-direction. 


As a result of relaticn (3.3.5), we have 
Ves) AVC) Re = MN), (3.3.6) 


Similar to relation (3595); we van obtain the following 
equalities 


SAU) Fe + Sal) FL =N [Se (t) con p+-5a (0) sin A] = NS, (1) 


Sit) Fe + SU FE — NSE (Uh eoew + Soe sin py] = NS (0). (303.7) 


In them angle u (Pig. 11) is determined by means of equations 


H = Neosp, FX = Nain, (3.3.8) 


and the constant factor N — by equality 
ae an 
wey (4) +B). (3.3.9) 


In right sides of equations (3.3.7) 5,08) — projection of the 
vector of the apparent path of the missile on the so-called hp- 
direction, which forms the constant angle p with axis € (Pig. 11), 
and Site) — current calculated value of this projection. By analogy 
with equations (1.6.7) and (1.4.7), 


t se 
SA =. \req@ae- homer. CS rduaet) 


Let us turn mow again to equation (3.2.13). Taking into account 
equalities (3.3.5), (3.3.6) and (3.3.7), we obtain 


MV.3 (0) + N [Sp (t) — Sp (t)] = MV’). (3.3.11) 
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Fig. 11. 


The right side of the last equation is a constant determined 
for each specific case of flight of the missile. This magnitude 
correct to the constant factor M is equal to the calculated value 
of projection on the A-direction of the vector of apparent velocity 
of the missile at the calculated instant of the switching off of its 
engine t = o*. If in equation (3.3.11) we introduce another designa- 


tion! 


Pa 
. mV ee (3.3.12) 


)’ 
2 ° 
(a5) + (se-) 
then it can be reduced to the following final form: 


Va (t) + PLS. (6) — 5, (2)] = Vi (6°). (3.3.13) 


For the construction aboard the missile of the left side of the 
ballistic equation (3.3.13), an additional integratior of current 
readings of only one integrator of accelerations is necessary with 
the axis of sensitivity oriented in the u-direction. Furthermore, 
the reproduction of the calculated values Sift) of this repeated 
integration is necessary. 


ee rE 


‘Let us note that coefficient p, similar to factor 1/1 in equa- 
tion (1.6.11), has a dimensionality opposite to time (in particular 


Soe. 
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§$ 4, Construction of the Ballistic Equation witn 
the Help cf a Special Meter of 
Apparent Accelerations of 
the Computer 


Let us examine another conversion of the ballistic equation 
(3.2.13), as a result of which the possibility of solution to the 
problem of inertial control of flight range of the ballistic missile 
hy metns of a single newtonmeter is explained. Its readings in the 
computer should be integrated after preliminary multiplication by 
a certain assigned time function. The axis of sensitivity of the 
newtonmeter in turn change its orientation of the nonrotating system 


vf coordinates according to the acsigned law. 


Let us turn, first of all, to the conversion of the first two 
corresponding terms of the ballistic equations (3.2.13). According 
to equations (3.2.4) and (3.2.8) and also the known Cauchy equation 
of the conversion of the repeated integral into a simple one contain- 
ing’ the upper threshold in the form <f the parameter in the sub- 


integrai expression, we have 
te ct ‘ . 
SO = (% (t)dt = (Gag (ey art =\« — vt) ag (x) de, (3.4.2) 
e ee 


where t — new variable of integration. Representing here the 
difference t - t in the form 


t—1t= (¢— 0%) + (0° — 2), (3.4.2) 


we cbtaln equality 
4 ¢ ; 
Se(t) = (t— 9°) { ae(t) dv + ( (0° — eh g(t) de. (3.4.3) 
e e 


heturning in the last integral of the right side of this 
equality to the initial variable of integration £ ani taking into 


account equation ©3.29.4), we obtain 


t 
Su) = (9M) + fo — Deine (3.4.4) 


Let us note further that correct to smallness cf the first order 
inclusively we can assume that 


. ors echt ; the 
u(t) = Ke (9) + (¢— 9) Ve (o). (3.4.5) i 
lee 
eon 
Equating to each other differences of the left and right sides 
of the last two equalities, we have 
t 
KOS) = Jie Nagnar— Sie) + ; 
+(¢— 0°) 1Me(t)— Valo). (3.4.6) 
equa 
The laSteterm of the right side of equality (3.4.6) has a ance 
second order of smallness and can be omitted. Actually, asa tion 
Mesaj. of tue proximity cf the real anu calculated instants 1 = ¢ ak 
ang t = of of the switching off of the engine quantity Veit) differs 
Little from Ve(o®), and this latter is distinguished from veto") by 
a small isochronal variaticn Svg(o#). Therefore the difference 
Ve(e) - ve(o®) is a 3mal] magnitude, which in equation (3.4.6) is 
Multiplied {r turn by anotner small difference t - a. 
Thus, correct to smallness of the second order 
: . 
Se (4) — St) = U0" — ty ag (t)dt — (3°) (3.4.7; 
. ® 
in perfect analogy 
‘ ‘ 
Su(t)— Sa(t) = \(3 — Nag (At — Se’). (a3 
$ . 
ang) 
Substituting differences (3.4.7) and (3.4.8) into equation ane 


(3.2,3%3) and replacing in it, furthermore, quantities v(t) and vice) 
Ly their reprecentsations, according to equations (3.2.4), we obtain 


after slapli fications 
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(+ —9 Fant [E+ —o Fen} = 
= M+ VO) + SUS) 1 FSi (0'). (3.4.9) 


The right side of the last equation for the specific flight of 
tne rocket is a certain constant, which we designate by the letter 

dn accordance with equations similar to (3.3.6) and (3.3.7), this 
constant can aisc be represented in the fourm 


Cm EVO) + eV) + HO) + FSO) = 
w MVi(3°) + NS, (9). “ (3.4.10) 


Conversion of the integrand expression of the left side of 
equation (3.4.9) can be produced by the same methods which were 
incorporated in the previous section. Let us present in this equa- 
tion alternating factors befcre projections of the apparent accelera- 


vibe a, il and a(t) in wre following manner; 


i 
’ 





RTO K(‘) 00s «(!), 


(2.4.11) 
& + OF K(t)sin x(t). 
| tn By corduiee with the iast equalitles, 
K(!) = V [+ e—o8] +[a, t+ -9F]. 
ny (3.4.12) 


The second equaticn (3.4.12) determines «xft) changing with time 


aligit: between the so-called «x-direction and axis &€ (Fig. 1?}. 
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Fig. 12. 


Using equations (3.3.3), (3.3.8) and (3.3.12), expression 
(3.4.12) for variable X(t), after comparatively simple conversions, 


can also be présented ‘in the form 
R(th=-M YT + 2p (oe —8) cos (up —¥) + ps — OF. (3.4.13) 


Let us produce in the left side of equation (2.4.9) replacement 
of coefficients according to equalities (3.4.10) and (3.4.11). we 
obtain the relation 


‘ . 
1X (6) [04 (t) cos (t) + a4 (t) sin x (t)] dt = C, (3.4.14) 
Z a, 
where in integrand expression the sum 
Ox (f) cos x (t) 4 ay (t) sin x (t) = a, (1) (3.4.15) 


is the projection of the apparent acceleration of the missile towaru 


the mentioned variable x-direction (Fig. 12), and constant C is 
As a result we obtain the ballistic 





expressed by equation (3.4.10). 
equation in the following form: 


t 
JE Mena =c. (3.0.26) 


(3.45.46) abvcard the 


| 
| 
sonstruction of the left side of equation 


missile :equires the presence of highiv succurate folicweup syetems 
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located on a gyroscopic stabilizer for the change in orientation of 
the axis of sensitivity of the meter of acceleration according to 
the assigned program. 


Of course, approximation methods of the formation of the left 
side of the ballistic equation (3.4.16) by means of the usual integra- 
tors of accelerations, for example, by means of replacement of func- 
tions Kft) and «(t) by their certain mean values. 


§ 5. Inertial Control of Range of the Ballistic 
Missile by Means of Longitudinal and 


Standard Integrators 
of Accelerations 


In systems of control of ballistic missiles one can use the 
so-called standard integrators of apparent accelerations, the axes 


of sensitivity of which are guided in plane En at the assigned angle 


to the axis € according to the program of flight. The control of 
Tuigi controts is producea in this case so that together with the 
fulfilling of the assigned program of the change in pitch of the 
missile (i.e., change in angle between its longitudinal axis and 
axis €) the reading on the standard integrator would be reduced to 
zers. The angle between the axis of sensitivity of the standard 
integrator and axis of the missile itself is selected close tou a 
stratght line. 


The use of the standard integrator of the apparent accelerations 
in the system of control of motion of the missile leads to stabiliza- 
tion in the assigned direction of the resuiting force of thrust of 
the engine and aerodynamic forces acting on the rocket. Thus, the 
direction of the vector of the apparent acceleration of the missile 
(hewe ver, of course, nol its magnitude) 1s stabilized. With the 
kiuown approximation, being distracted from errors of the system of 
control, it can be considered that the standard integrator controls 
the flight of the missile so that the projection of the apparent 
nuceeleration on the axis of sensitivity of this integrator would be 


equal to zere. 


? 





Sok 


Let us present in the ballistic equation (3.4.16) of the 
previous section quantity a, ft) in the form of the product 


a. (t) = a(t) cos z(t). (3.5.1) 


Here y(t) is the angle which forms the vector of the apparent 
acveleration a(t) with x-direction inclined at angle «xft) to axis 
& (Fig. 13). 





Fig. i3. 


Function «(t) is determined by the second equality (3.4.12), 
and the direction of the apparent acceleration is changed according 
to the assigned law. Because of this the magnitude of angle y(t) 
snoulid be considered as a known function of time. 


Producing in equation (3.4.16) the replacement of quantity 
a,(t), according to equation (3.5.1) we will arrive at the follow- 
tng modification of the ballistic equation: 


Somemaeac, (3.5.2) 


in which the aiternating coefficient 


Q(t) = A(t) cos 7 (8) (3.5.3) 











should be considered the known time function for the given cal- 
culated case of flight of the missile. 


For construction of current values of the left side of equation 
(3.5.2) aboard the missile just as earlier (see § 4 of this chapter), 
integration of current values of the apparent acceleration, multiplied 
by the assigned time function is required. However, now there is no 
need in special high-precision equipment for the change in oriertation 
of the axis of sensitivity of the meter of the apparent acceleration 
of the missile. Actually, if this newtonmeter is rigidly connected 
to the standard integrator so that axes of their sensitivity vould 
pe perpendicular to each other, then with the proper accuracy of the 
control of flight of the missile the apparent acceleration of the 
latter (more aecurately, apparent acceleration of location of these 
instruments) wiii te wholly directed along the axis of sensitivity 
of tne mentioned newtonmeter. 


“4th some Jess of accuracy of flight of the missile the rnewton- 
meter can be fastened directly aboard the missile, the guiding axis 
of its sensitivity in parallel to the longitudinal axis of the 
miscile. The iatter 1s ueflected from the direction of the apparent 
acceleration of the missile, i.e., from the direction of the resulting 
force of thrust of the engine and forces of aerodynamic actions, as 
a rule, at small angiec. Because of this, approximately to lay 


eo (t) 0 (8), (3.5.4) 


where a(t. ~ projection of apparent acceleration of the missile on 
its longitudinal axis. 


lf,in accordance to the last approximate equality, we produce 
the appropriate replacement in equation (3.5.2) and, furthermore, 
subctitute function Q(t) by its certain mean value q, then we will 
obtain the approximate equation 


8 
QS an(edt = QVC (3.5.5) 
e 





‘for the determination of the instant of switching cff of the engine 


of the missile. Here expression 


% t ~ 
Vos) = S cateyat (3.5.6) 


is the current reading of the integrator of acceleration with the 

axis of sensitivity parallel to the longitudinal axis of the missile.! 
Such an integrator is called longitudinal. As is known, it was used 
in the guidance system of the German missile V-2. 


Let us return again to equation (3.5.2) and examine additionally 
certain approximate methods of construction of its left side on the 
basis of current readings 


t 
Minn Jewma (3.5.7) 


of the integrator of accelerations with the axi: of sensitivity 
located perpendicular to the axis of the standard integrator. Con- 
sidering the latter equality the left side of equation (3.5.2: is 
converted to the form 


e ‘- ae 
yecnacnarm (Qh a (3.5.8) 
e 


and can be integrated by parts. As a result, taking intc considera- 
tion that V(O) = 0, we obtain, according to relations (3.5.0) and 
(3.5.2), the ballistic equation in the following form: 


$ . 
@mron--ferw Pier: ¢. (3.5.9) 


‘Let us note that quantity v(t) in equation (3.5.¢) because 
of the variability of direction of the longitudinal axis is not the 
urojection on this axis of the apparent velocity of the missile in 
‘he nonrotating system of coordinates (for more detail] on this ser 
§ 4 cf Chapter I). 


Nan 


The construction aboard the missile of the left side of the last 
equation can already be produced by means of only one integrator of 
acceleration. 


Let us note in conclusion that the same method of integration 
by parts can be used for constructing the left side of the ballistic 
equaticn (3.4,16). 


‘+ camming Mala ile iA na pts ems = 





mem 
By 
be 
the 
and 
CHAPTER IV 
INERTIAL CONTROL OF THE RANGE OF PLIGH! CF THE oe 
BALLISTIC MISSILE ALLOWING FOR THE 
CHANGE IN THE FORCE OF GRAVITY tie 
€xp 
§ 1. One Method of the Solution of Dirferential whit 
quations by Which Isochronal Variations of ee 
Coordinates of the Missile with Its 
~~WMotion in the Powered-Flight of 
ection Are termined oe 
eng 
in the previous cray-ter the basic ballistic e;uation (2.7.6) was Tak 
transformed to such a rorm at which its lert side cuuld be constructed iso 
abcard the missile with the help of integrators of accelerations. In we 
this case iscchronal variations of coordinates and projections of 
velccity of the missile, which enter into the composition of equa- 
tion (2.7.6), were expressed by the approximate equations (3.2.6) 
and (3.2.2). This corresponded to the neglect in differential equa-~ 
tions (3.1.3) of terms containing as factors the isochronai varla- 
tions S&(t) and é6n(t) tnemselves. As was already indicated in cl:cpter 
III, the mentioned terms of equations (3.1.3) take into account dur- 
iKg@ the determination of current coordinates cf the rocket changes in am 
acc@leration of the force of gravity because of the noncoincidence of 
its actual motion with the calculated. Error appearing from such a 
simplification of equations (3.1.3) in the determining of the range 
of flAght of the missile is small only at rather small deviations in che 
its real motion from the calculated. Otherwise it 1s necessary to of 
solve approximately aboard the missile differential equations (3.1.3) 
3llowing for terms reflecting the effect of the change in accelera- 
tior. af tne force or gravity. sol 











Below for the solution of equations (3.1.3) a similar approxima- 
tion method is used, and it was shown to the author by corresponding 
member of the Academy of Sciences of the Ukrainian SSR Yu. D. Sokolov.? 

' By means of this method the problem of inertial control of range can 
be guite accurately solved by the means given in Chapter III with 
the help of the application of standard integrators of accelerations 
and the simplest ccmputers. 


The idea of the method consists in the conversion of equations 
(3.1.5; to their equivalent integral-differentia. and integral rela-~ 
tions with sursequent replacement under the sign of integrai of 
expressions of jesired functions by their simplest approximations, 
which satisfy, however, initial and final conditions of the problem. 
In accordance with this, let us integrate the right and left sides 
of equations (3.1.3) with respect to time from the initial instant 
t = O up to the unknown until the instant of switching off of the 
engine t = 6, which provides the hitting of the missile on the target. 
Taking inte ace out in this case initial covditicns (3.1.5) for 


Sn a, 


cyan eli alee al cabs te 25 


isochronal variations or projections of the velocity of the missile, 
we obtain two integral-differen:ial relations 


e a 
FE) os bug (a) = O¥e(3) + Vaart Fone 


a, 
dq Ontedt, (3.1.1) 





e . e 
(a st 
#319) = bu, (3) 2d 4(3) t\e bdr 


in which, similiar to equalities (3.2.2), 
We(3) = [Bag at ov.) = §teg¢nat : (4.1.2) 


the isvenuronal variations of projections of the apparent velocity 


of the missile on axes & and n at instant t = o. 


‘See, for example, Yu. D. Sokolov. One method of observed 
solved linear integral differertial equations. — Dop. AN URSRF, 
1955, No. 2. ‘ 
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If, however, the right and lef* sides of equations (3.1.3) are 
integrated with respect to time twice: one time within the limits 
of t = 0 to the current instant t and the second — also from t = 0 
to the instant of switching off of the engine t = o, then, taking 
account initial conditions (3.1.4) and (3.1.5), we obtain an addi- 
tional two integral relation: 


= 


e. 
N {se onion, 


on 866 


8E (6) == 5Sz (5) + 


re 


fe 
a Rae + 


es oi 
8n(a) = 884 (0) + We BE (t) de + \\ ae anne, (4.1.3) 


where in turn, in accordance with equalities (3.2.6), 


BS_ (0) = Gareqnar - § ay (2) d¢*, 


(4.2.4) 


¢ 
8S, (a) = Goren = ( (da, (8) dt? 


{isochronal variations of projections of the apparent path cf the 
rocket, which refer to the instant t = o. 


In accordance with the mentioned method, one should substitute 
in the right sides of relations (4.1.1) and (4.1.3) functions 6&(t, 
and 6n(t) by their approximate representations, which turn into zero 
at the initial instant of time and respectively into s&(o) and én(o° 
when t = o. Then values of quantities 6&{0', in(o), du,/o) and 
bu (a), found in accordance with the mentioned relations, prove to 
be, as a rule, more accurate than those calculated according tc 
equations (3.2.2) and (3.2.6), founded upon approximate equalities 
(3.2.1). 


it is possible to take as approximate representations 6€(-:/ and 
fncett in right sides of relations (4.1.1) and (4.1.3), for exampie, 
tly i: .inear approximation 








tak OS, n= anise, (4.1.5) 


where o — instant of time for which the determination of desired 
functions is produced; in this case o — instant of termination of 
the powered-flight section of flight of the missile. 


The best results should be expected with the quadratic approxima- 


tion of the form 
BEC) = BEA)-Z. —-dn(t) = Bn(0)-4. (4.1.6) 


The reason for this is that functions (4.1.6), unlike functions 
(4.1.5), satisfy simultaneously initial conditions (3.1.4) and 
(3.1.5), which concern both variations themselves of the coordinates 
and their time derivatives at the instant t = 0. Functions (4.1.5) 
do not satisfy initial conditions (3.1.5). 


Leading to am even greater accuracy should be the assignment 
of functions 6€(t) and 6nft) in the form of the following polynomials 


of the third power 


8E() m= 8E (0) (3-H — 2G) + dme(s)-0 (G — F): 


dn (t) = dn (0) (3-4 —2 5) + du (0)-0(S — 5). (4.1.7) 


These polynomials not only satisfy initial conditions (3.1.4) and 


(3.1.5), but, furthermore, when ¢ = o turn, respectively, into 5&(o) | 


and Sn(o) and their derivatives — into Su, (3) and Su (o). 


& 2. Approximate Solution of Differential Equations 
for Variations of Coordinates Using the Quadratic 


Approximation and Also the Approximation in the 
Form of Polynomials of the Third Power 


Let us examine from the beginning the solution of differential 
equations (3.1.3) which corresponds to the quadratic approximation 
ci.4d.6). Suverituting into right sides of relations (4.1.3) 


6b 
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expressions for S&(t) and é6n(t), according to equations (4.1.6), 
we obtain 


" 8B(o) = 8S (0) + goth He tis +t Oe aan, 
; ee 


of, , os : 
8n (0) = IS, (0) +P ae ae + UN Se ae: (4.2.1) 
‘ e e . 
Let us introduce here designations 


et ot 
PaO) = ae \ Ed: an 0) = ar GE ae 

“ne of (4.2.2) 
Bao) = or | a ee bento) = ae) ee, 
x 


where, because of relations (3.1.6), 


Ins (0) = ha (0). (4.2.3) 


As a result let us arrive at the two algebraic equations: 


88 (6) = &Se (6) + Aen (0) 8E(0) + Aes () 8n(0), 
8m (0) = 854(0) + hae (0) 88 (0) + hee (0) 8n (0) 


with respect tc the desired quantities 6&(o0) and én(o). 


Values cf coefficients heels), he, (0) = hye (o) and hin lo) with 
an accuracy sufficient for practice can be taken with argument o = 
2 of, Actually, the real duration of the powered-flight section o 
only by a small magnitude is distinguished from the calculated o*. 
Isoenronal variations 6&(0) and é6n(o) should also be considered as 
small magnitudes. Consequently, the error which occurs from the 
ment toned replacement in coefficients h,,(o), he, (0) and hoy) s 
which stand in the right side of equations (4.2.4), of the value of 
arrumnent o by o#, has the second order of smallness. 


69 








Let us designate coefficients (4.2.2) when o = of by the same 
letters but without Indication of the argument. With the made 
simplifications and replacement of the designation of argument oa by 
t equations (4.2.4) take the following form: 


OB (6) mm O5g (8) + AqydE (t) + dead (2), (4.2.5) 
my (t) am Sy (2) —t Age BB (t) + Auydn (t), 


where ligge hen = Ane hen are quantities known for the given cal- 
culated fiight of the missile, which are small in comparison with 

unity (see § 3 of this chapter). Because of this the totality of ! 
equations (4.2.5) is must convenient to solve by means of iterations, 
taking 6€(t) and é6n(t) in the right side of equations (4.2.5) the 


4 
basic expressions represented by equations (3.2.6), namely ' 


BE (t) we OSe(¢), Om (t) = 5, (6). (4.2.6) 


by substituting these expresctons intc right sides of equations 
(4.2.5), we arrive at the following equations for current isochronal 
variations sf coordinates I8(t) and é6n(t): 


E(t) mm (1 + gy) d5g (6) + AgadSs (¢), 
my (8) we hag dSe (1) + (1 + Ang) 5, (0). eee 


it is easy tc be convinced that subsequent iterations of equa- 
tions (4,2.5) ieai accurately to the same equations for €&(t) and 
én(t), if only in the calculations we drop terms with squares and 


products 4 h 


Es, “Eq Nee SNe Man: 


To search for variations of projections of the velocity of the 
misslie Cup lt) and Sui(t), let us substitute into right sides of 
relations (4.1.1) expressions for 6€(t) and énft), represented by 
equations of the same quadratic approximation (4.1.6). We obtain i 


for the desirea magnitudes these equations 








so 


(4.2.8) 


i= Sls 


Sue (o) = SV (0) + SE (0) wy 3g Pdt + bn (0) a) = titdt, 
Sa 
#\e dt. 


du, (0) = oV,(0) + worse Se Sp dt + bn(o) 
Let us introduce designations 


e 
fu (0) = 5 mye + fta(6) = ~ s)he, 


(4.2.9; 


sal) = lag Mat seal) = See a 


Then equations (4.2.8) are presented in the form 


“Bag (3) = 8V'g (0) + sez (9) 88 (5) + Fen (2) 5n (0), 


(4.2.10) 
8u, (0) = 8V,(0)+ gaz (5) 5E(5) + gan (0) 8n (6). 


Similar to the previous in equations (4.2.10), instead of 
coefficients Bee lo), Ep, (0) = Enel?) and By 6)» it is possible to 
take thelr values at the calculated instant of the termination of the 
powered-f'light section of flight of the missile o*. We designate 
these values, respectively, by Begs Ben = Bne and Ey, Let us sub- 
stitute further into equations (4.2.10) expressions ine -5) for 
variations of coordinates 6&€(ca) and 6n(o) and disregard the products 
of coefficients g and hk with arbitrary indices. Ultimately we 
arrive at equations for the desired isochronal variations of velocity 
of the missile up (a) and du, (a) with*respect tc immobile system of 
coordinates ng. Changing in them the designation of argument o by 
L, we obtain 


Sag (t) me SV g(t) + Kee bSe (t) + 8e.555 (2), 


(4.2, 
Bang (t) = BV a(t) + 2ed5e (4) + KardSs (2). ay 


Equations (4.¢.7) and (4.2.11) are the appru«imute solution of the 
totality of differential equations (3.1.3) in a “orm similar for the 


fed, 





application in systems of inertial control of flight range of 
ballistic missiles. 


Somewhat more bulky are equations for isochronal variations of 
coordinates and projections of velocity of the missile, if we take 
polynomials (4.1.7) for the initial approximations in equations 
(4.1.1) and (4.1.3). As a result, similar to the case of the. 
quadratic approximation, it is possible to arrive at the following 
-equations for variations of coordinates: 


88 (t) = ({ + Dre — 2fex) Se (t) + (Bhea — 2fes) 5S, (t) + 
+0" (Fax — hagy) BV (1) + 5° (Jen — ee) OV (8), 


3m (8) = (Veag — Daz) OSe (0) + (1 + Bega — Ban) 85e (t) + 
$5" Ua, — eae) Ve (t) +9" (/aa— hae) V4 (t), (4.2.12) 


and for variations of projections of the velocity, correspondingly 


Bae (t) = (Beez — Qige) BSe (¢) + (Been — Bea) BS, (8) + 
+11 + 9° (ieg— gee)] OMe () + 9° (ing — 8:5) OV (0), 
Beg (t) = (Bean — 2hne) BSe () +- (BB aq — Vhan) OS, () + 
+ 8° (Eg — Gra) OV (0) H+ [1 A 5° Finn — Bae)} O00 (0). (4.2.13) 


In equations (4.2.12) and (4.2.13), besides designations already 
introduced in this section, quantities Urps tens tee Try @ and Tepe 


Jen? Ines dnn are values of integrals 


e © ' 
, we a: edt, i (0) = 4 \ 52 ear, : 
bye (9) a) » fas 85 an 
(4.2.14) 2 
ign (0) = tna (0) =e +\aea 
ured 7 et ‘3 ‘ 
pine EL Mera, iat = deft Bee, 
ee 
(4.2.15) 
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hale) = fa (9) = Hay al \ 





which they take with argument o equal to «#, i.e., at the calculated 
instant of switching off of the engine. 


§ 3. Simplification of E uations for Isochronal 
Variations of Coordinates and Projections of 
elocity of the Miesile at Extensions of 


the Fowered-Pilght Smaii in Comparison 
with the Radius of the Bart 


Equations (4.2.7) and (4.2.10) and also (4.2.12) and (4.2.13), 
obtained in the previous section for variations of coordinates and 
projections of velocity of the-missile, can without great harm for 
their accuracy be represented in simpler form, which dces not require 
preliminary calculation of integrals of the type (4.2.2), (4.2.9), 
(4.2.14) and (4.2.15). Actually, partiai derivatives of the accelera- 


tions of the force of the earth's gravity zg az zt taking 


part in these integrals, are functions of moving coordinates of the 
calculated motion of the missile and, consequently, finally — time. 
However, with a smali extension of the powered-flight sectior -' une 
trajectory of the missile, in comparison with the radius of the 
2arth, these functions can be substituted by their values, which 
refer to the point of the missile launch, and the derivatives them- 
selves can be calculated on the assumption that the earth is a sphere 
with a radial distribution of density. 


The last assumption, if one considers equations (2.5.4), leads 
to the following expressions for the mentioned derivatives of 
accelerations of the force of gravity: 


R--F-3$).  H--Fl'-3H), 


jean te. (4.3.1) 


Here &, n and ¢ — coordinates of the arbitrary point of space 
and p — its distance to the center of the earth. 








For the point of thé missile launch (see § 1 of Chapter II) we 
have 


B—=F=0 and wpa, (4.3.2) 


and, consequently, at this point, according to equations (4.3.1), 


he o/, o 
A Be BWo 
$e = daw, i (4.3.3) 


where there iz: introduced designation 


ve VF. (4.3.4) 


In the theory of gyroscopes and inertial navigation quantity 
v, determinec by equation (4.3.4), is known by the name of Schuier 
frequency (v = ¢.90125 l/s). 


Ysing the approximate expressions (4.3.3) for derivatives 
ey @ al Pe : 
By BHF and st ir equalities (4.2.2), (4.2.9), (4.2.14) and 


(4.2.15), we arrive at tie following simple equaticns: 


her (8) = — FF. ee (0) = hea) = 0, Aenlo) = HF, 
cu) = — FF. tenlo)— 5a), trio) = 7, (4.3.5) 


and als¢c 
ha) = — SF. ke) fal) = 9, info) =e, 


ig()= — SE, tn) ia) = 0, inne ZF. (4.3.6) 


Having assumed in these equations o = o#, one can present expressions 
(4.2.7) and (4.2.11) for 4tsochronal variations of coordinates and 
projections of velocity of the missile with the initial quadratic 
approximation (4.1.6) in the following final form: 
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a) =[1— “SF asco, 
én) =[1 + 22F Jasscn, 
Sug (t) = OV e (t) — SF 85400), (4.3.7) 


butg (t) = 8V g(t) + = 85, (1). 





Similarly, for the case of the. initial approximation of functions 
6E(+) and Sn(t) in the form of polynomials (4.1.7), on the basis of 
“equations (4.2.12) and (4.2.13), we obtain 


aio = [1 — AEF as. + AEF oven, 
b(t) = [1 + de 6] 85.0 — 240)" 47, (0), 
ug (4) = [1+ “EEF ] 0% 1) — F350, 


du, (f) = [1— AEE] OV a(t) + v0" 854 (0). (4.3.8) 


Equations (4.3.7) or (4.3.8), wnich approximately consider the 
change in action on the missile of tne force of the earth's gravity 
‘necause of the tmotion of the missiic not according to the assigned 
ealculated law), should be more accurate than similer equations 
(3.2.2) and (3.2.6) given in Chapter III. 


§ 4, Construction of Ballistic Equations in Which 


the Effect of Changes in _Eitect of Chenkes in Accelerstion of the 
Force of the Earth's Gravity i 
Approximately Taker, Account 


In the previous paragraph equations (4.3./) and «:s0 (4.2,8) 
were obtained for the determination of variations of seoordinate- and 
projections of the velocity of the missile with °~ ippr: ximate 
calculation of the effect of the change in the f-rce of the earth's 
gravity with motion according to the law which ts somewhat distin- 
guished from the calculated. The mentioned equations allow nore 
accurately, in comparison with expressions (3.3.13) and (3.4.16; 


1S 
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obtained in Chapter III, constructing the left side of the basic 
ballistic equation (2.7.6) aboard the missile, using readings of 
the same integrators of accelerations. For this we will use in the 
beginning equations (4.3.7) and substitute into equation (2.7.6) 
expressions for é6&(t), én(t), Su, (t) and’ 6u,(t) determined by them; 
we obtain 


850) {ft — | e- > alt 
aS at 


+asn(o {ft + 2EF] A 4 Ae Ah 
FOR THOR, = 
~—¢— 39a Oa tee tel: (4.4.2) 


Further, similar to that which was done in § 2 of Chapter III, 
jet us present variations of projections of velocity $V, (t) and 
SV_ft) in the form of differences (3.2.3) and use equations (3.2.12). 
Ir tnfies tuse gust as in the conversion of equation (3.2.10), terms 
standing in the right side of equality (4.4.1) are reduced with 
winilar terms of its ieft side. Taking into account still equations 
'g.¢.74, we arrive at the equation 


1Se (8) —Se(#)) Ae + 154(— SAA + Ve OF + 
+ Vathgge = VES) Sot Va (44.2) 


in which 
Mere a vet 
"Ot EE OE ST” Guy * 
VSP at | We" al 
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(4.4,3) 
A, = La tie 


Equation (4.4.2) differs from equation (3.2.13) only by 
‘oe fficlents of ftsochronal variations of projections of the apparent 
path So, (t) and 65, (t). Such pee cients in equation (3.2.13) are, 
respectively, p day derivatives = and = and in equation (4.4.2); — 


+ 


APA i UR i it spill a edt 


ite 


are quantities Ay and Ay distinguished little from them. Therefore, 
subsequent conversions of equat on (4.4.2) can be produced exactly 

us in § 3 of the previous chapter. Then instead of equation (3.3.13), 
we will arrive at the following: , 


Va (8) + PIS, ()— Sp (¢)] = Vi CS). (4.4. u 


The left sides of equations (4.4.2) and (4.4.4) are constructed 
aboard the missile with the heip vf the same technical means as for 
equations (3.2.13) and (3.3.13), i.e., two integrators of accelera- 
tions and the appropriate computer, which includes elements of 
repeated integration. In this case the axis of sensitivity of one 
of the integrators is at the same constant angle A to axis € as ir 
the construction aboard the missile of the left side of equation 


(3.3.13). As regards the axis of sensitivity of the other integrator, | 


the latter should be inclined toward the axis & at an angle of u', 
which is little distinguished from the corresponding angle yp. Similar 
to equality (3.3.8) the magnitude of angie u' is found from relations 


Ag= N'cosp’, = Aye N’sinp’, (44,6) 
where, of course, 


W’ = V (Ag) + (Aa). (4.4.6) 


In turn the coefficient p’ is distinguished little from coefficient 
_p of equation (3.3.13) and is determined by equation 


va >. (4.4.7) 


Here quantity M is the same as that in equation (3.3.4). 


Possible also is the conversion cf the ballistic equation 
(4.4.4) to the form 


fr wamaac (4.4.8) 
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similar to the form of equation (3.4.16). The use of equation 
(4.4.8) allows solving the problem ebout inertial control of flight 
range of ballistic missiles just as in § 4 of Chapter I1I by means 
of cone meter of acceleration with an alternating direction of the 
axis of sensitivity. In equation (4.4.8) the variable coefficient 
K'(ti is formed t; the same equations (3.4.11) as coerficient K(t) 
in equation (3.4.16) but with the replacement of angle «(t) by angle 
Kk' ft) and darivatives x and z » respectively, by quantities A 


and A‘,, 
n 


The use in the transformation of the basic ballistic equation ! 
(2.7.6) of more compliex equations (4.3.8) for coordinates and pro- 
jecticns of velocity of the missile leads, naturally, to more bulky 
results in cenparison with equations given in the beginning of this 
section. Instead of equation (4.4.1), in this case we arrive at 


the more complex equation 


wie REY OO obs Rh ante setntenastes oo ta 


as (25 Es 
+5.{ft + aver] tye F} + 
rel BE afte E+ 
ore( eras 2Fr)2) 


=—(t-o) [ROR + eS. (4.4.9) 
in turn equation (4.4.2) 18 replaced by the following: 
1Su (6) — SUC AG + 156 (0) — SA) AN + HE) BU 


+.) Bi — Vito (B— 2) — VBL — ZS) = 
—Mee) K+ Vale) R.- (4.4.10) 


In it these designations are introduced 


A [1-5] & - FH 











Avm [t+ ES 4 0 2, 
in ett ]z. (4.4.11) 
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Similarly, equation: (4.4.4), when using equations (4.3.8), is 
replaced by the following: 


Vy (t) +p” [Spe (t)— Spe (t)] = Vie (0°) + 


t—<s° 


+492 (Be— gr.) e000) + ( (8.—Z)ae]. (4.4.12) 


Here A" and uw" are directions of axes of sensitivity of two integra- 


and u-directions, introduced in § 3 of Chapter IIl. Angles \" and 
u", which these directions from with axis £, are determined by 
relation 


Ag= N'cosp, A= N' sing, 
By = Moos), By= M sind. (4.4.135 


Coefficient p"” in equation (4.4.12) is expressed by equation 


pase (At (4a? “ok q 
mY) wre eee 


the construction of the right side of the ballistic equation in 
the form (4.4.12) aboard the missile is complicated, as compared to 
the case of the equation of the form (4.4.4), by the necessity of 
introduction into the computer of an additional term iinearly 
jenendent on time. Finally, subsequent conversions of equation 
(4.4.12), specifically, to the form similar to the ballistic equa- 
tion '4,4,8) and also (3.5.2) are allowed. 


'n ecnciusion let us note that the selection for the specific 
aveotem of control of the flight range of one of the baliistic 


tors of accele ations, which, respectively, differ little from \- 
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equations given in this and previous chapters should be produced on 
the basis of estimates of quantities of the so-called systematic 
errors inherent to the system of inertial control of the given type 


of rocket. Systematic errors are understood usually as those errors 
in the determining of flight range of the missile, which appear 
exclusively because of simpiifications made in the derivation of 

the appropriate ballistic equation. Because of this, even with the 
accurate switching off of the engine of the rocket at the instant 
indicated by the ballistic equation and with the accurate operation 
of al] remaining instruments and devices of the system of inertial 
control, the actual flight range of the missile can be somewhat dis- 
tinguished from that assigned according to the calculation. As 
regards to the technical difficulties of the construction of systems 
of inertial control of flight range of ballistic missiles, they are 
included in the manufacture of meters and integrators of the apparent 
accelerations (1.e., newtonmeters and impulse meters) with extremely 
small instrumental errors and gyroscopic instruments with minimum 
servicing of their stabilized axes relative to directions on fixed 
star. ‘The same refers to the immediate switching off of the engine 
arn the signal of the achievement by the ballistic function of a value 
2orresponding to the assigned flight range of the missile. 
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CHAPTER V 


CONTROL OF LATERAL MOTION OF THE BALLISTIC MISSILE 


§ 1. Lateral Deviation of the Missile from the Target, 
Expressed in the Form of a Function of Changes in 
Parameters of the End of the Powered-Flight 
Section in the Starting system 


of Coordinates 


in the solution of the problem about the elimination of lateral 
deviation in the missile from the assigned target, there appear 
additional aifficulties connected with the fact that cf the duration 
of the free-flight section of real motion of the missile, as a 
rule, 1s somewhat distinguished from !*s calculated value. Tn view 
of this, control according to the earlier assigned law of the motion 
of the missile in a direction perpendicular te the programmed pane, 
i.e., plane zy of the starting system of coordinates xy2 oor (in: 
another possible variant) to plane &n cf the nonrotating sy-ctem 
énc, in general does not provide the absence of lateral deviation 
in the missile from the target. Really, because of the laws of 
mechanics the missile on the createi part of the free-flight section 
of its flight is moved with minute deviations from a certain plane, 
which does not change its orientation relative to directions at 
fixed stars (the mentioned deviations are connected hasically wit? 
the nonsphericity of the earth). Therefore, if the duration of th: 
free-flight section is not equal to its calculated value, ther. 
because of the rotation of the earth lateral deviation in the 
missile fron the target’ continuously appears. However, with 





‘An excepticn ‘s the case of flight of a missile in the plane of 
tle equator ané from one pole on another. 
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small distinction in the actual motion of the missile from cai- 
culated this deviation should also be small. Correspondingly the 
control of motion of the missile in the powered-flight section in 
a lateral direction can be reduced only to the requirement of the 
fulfillment of an earlier assigned (calculated) law of motion in 
the projection on the z axis (or in another variant - on axis ©). 
However, for a more accurate elimination of lateral deviation in the 
missile from the target it is necessary to take into account how 
at current instant of time its real motion in a longitudinal 
iireectior is distinguished from the calculated. For this pvurpose 
one can use current readings of instriments of the control of 
flight range sf the missile. 


It is obvious that an especially accurate control of lateral 
motion of the missile is necessary only during a short interval of 
time when cutoff of its engine can occur in flight at the assigned 
range. At the remaining time of the powered-flight section the 
accuracy of control of lateral motion can be somewhat less. The 
duration of the mentioned interval of allowed cutoffs of the engine 
depend on the quality of control of motion of the missile according 


to the assigned program. 


Tn the starting system of coordinates xyz (see § 1 of the 
second chapter, Fig. 9) quantity b - of lateral deviation in the 
missile from the target, just as the range Tr 1s a function of 
coordinates x, y, 2 and projections V9 V,s Y, of its velocity 
relative to this system of coordinates at the instant of termination 
of the powered-flight section of the flight. Thus, 


Ce (5.1.1) 


The astronomical time, which corresponds to instant of termina- 
tion of the powered-flight section, does not play any role here as 
in equation (2.1.2), unless one does not take into account the 
negligible effect of such factors as the mutual arrangement of the 
earth and moon or the earth and sun. 
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With the replacement of arguments of function (5.1.1) by their 
calculated values r#, y¥*, 2*, v¥, v¥®, and Ds it should turn into 
-zero, since the lateral deviation in this case is absent. Thus, 


b(z*, y*, 2°, of, vy, vt) = 0 (5,722) 
and, consequently, 
Ab = b(z, Mr 2, Oxy Dy v,)—b(z*, y’, 2°, vw, vy, v) = 5, ( 5 . 1 . 3) 


f.e., the error in th2 lateral deviation Ab is eaual to the latera? 
deviation b. 


In accordance with the Taylor series for the function of manv 
variables, we have, correct to smalls of the first order inclusively, 


the following expansion: 


bm (z— 2) Ft y— YF — B+ 
+0, — 02) + (0,— 04) ee + (0-07 a (5.).h) 


Here derivatives =. 2 a oe we. Sie -. ure *aken at cal- 


or’ Oy’ Os’ dry? Ory OP, 
culated values of their arguments and, consequentiv, for each 
ealculated case of flight of the missiic are earlicr Kucwn quantities. 


‘They can be determined in a way similar t> the ballisti: toeffietents 
given in § 1 of the second chapter. 


’ 2, Lateral Deviation as a Function of Pa:meter: 
of the End of the Powered-Flight Section in 
the Nonrotating System of Coordinates 
In examining the problems of inertial controi ** °°\".era} 
motior of the ballistic missile, the starting system 7° ccordi- 1ves 
1s inconventent. Therefore, let us introduce the same worn riat one 


_ 


system of coordinates Ent as that in the examining of the the: |: 
'nertial control of range (see § 2 of the second chapter). Current 
coordinates x(t), y(t), 2zft) and projections v(t), vy ft) and v(t) 
“of tune velocity of the missile relative to the starting system of 


33 








coordinates xyz are expressed in this case by current coordinates 
E(t), nt), tft) ana projections up(t), ur (t) and ur (t} of the 
velocity of the misslle in the nonnobatene system ae coordinates 

Eng by equations (2.2.5) and (2.3.5). Entering into the mentioned 
equations also in evident form (through angle ¢ = Ut) is time t, which 
passed at the instant of launch of the missile. In view of this, 

in accordance with equation (5.1.1), the magnitude of lateral 
deviation in the missile from the target can be represented in the 


form of a function! 
b aw B(B, 1, 0, ug, Uy, Ue; 3) (5.2.73 
or ecordinates €, n, % and projections Ups Uns 


of the missile relative to the system of coordinates Enz at the 
instant t = o of términation of the powered-flight section and also 


uy of the velocity 


tvs very duration co. 


Calculated values of coordinates x*, y*, 2* and projections 
ve, ve, oe cf velocity of the missile relative to the starting 
ayanen tyz correspond to certain calculated values of coordinates 
E*, n*, ¢* and projections ues us, ut of velocity of the missile 
in the nonrotating system Ent. It is obvious that function (5.2.1), 
which represents the lateral deviation of the missile from the 
target, will turn into zero at these values of their arguments if, 
furthermore, we assume the duration of the powered-flight section 


is also equal to its calculated value o*. Thus, 
b(E*, n°, 0, we, us, Ue; O°) a = 0, (5.2.2) 


Expanding function (5.2.1) in Taylor series and being limited 
in the latter only by terms of the first order of smallness, we 
obtain, taking Into account equality (5.2.2), the equation 


PH E-— NSH — MZ +E S4 


FOOD Fe Mae HH — MZ +0 %, (5.9.3) 
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tere, fust us in § 4 of the second chartor, jeciznattons 
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are aotented. 
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in which partial derivatives a 2, $. me he and = are a 
constant coefficients, since their arguments &, n, &%, Urs Uns Uy ita 
and 9g are taken equal to their calculated values, i.e., respectively Sat 
E*, n*, c®, urs uss ue and a*, Magnitudes of these coefficients 
depend finally on the selection of a certain calevoated case of 
flight of the missile. In a way similar to the appropriate partial fune 
derivative of function 2, which represents the flight range of the sign 
missile, they are expressed by equations of the type ‘c.4.5) and This 
(2.4.6) in terms of partial derivatives =. =. =: ae me x turn 
and also quantities x#, y#, 3%, v#, vi, vi. TAS" 
a eee deri 
Equations (5.2.3), correct to squares and products of differences 
R—~B n—n C—O, ue — ue mh, mem o—o* , determines the lateral 
deviation in the missile from the target when magnitudes of coordi- 
nates and projections of velocity of the missile at the end of the 
powered-flight section are not equal to their calculated values. 
of la 
With real motion of the missile neither the instant «af switching of f 
off of the engine o nor its coordinates &, n and projection of makim 
velocity Ups ys which correspond tec this instant of time, are B.(t) 
Known earlier. Therefore, for elimination with the above-mentioned 
accuracy of lateral deviation in the missile from the target, a 
continuous change is necessary in quantities fit) and uy fe. charac- 
terizing it~ motion in a lateral direction during a certain time 
interval of the end of the calculated powered-flight section. The 
purpose of such a control is the continuous reduction to zero of the 
funetton missi 
nonro 
B= RO — ENE + 9 WF + V0 S + 
+ le ()— Wel ge + MO ‘ ; 
r + lee) ual et NS (S.7,u ecnsi 
I and x 


in the whore .:terval of possible instants of the switching -ff .f 


' the engine for a missile of the giver. system. Actually, functior 
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@(t), according to equations (2.4.1) and (5.2.2). turns into the 
magnitude of lateral deviation in the missile b if we assume in 
it argument t equal to the actual value of the instant of the 
switching off of the engine o. 


The current value of function @(t), which we will call the 
funetton of lateral deviation, can be used as the controlling 
signal in the system of control of lateral motion cf the missile. 
This system should operate in such a way that function 6(t' wouid 
turn into zero because of the corresponding change in current 
magnitudes of the coordinate of the missile cit) and its time 
derivative up(t). 


§ 3. Construction of Corresponding Terms of the Fune clon 


of Lateral Deviation Dependent Onl 
the Basic Motion of the Hisstic 


As was shown in the previous section, for the inertial control 
of lateral motion it is necessary to reproduce the curre:t value 
of funetion 8(t) aboard the missile. Tne totality cf tems (5.2.4) 
making up this function is decomposed into two groups BL t) and 
6 .'t). The first of them is the sum 


Aw =) — PS umes iz + 
aa a 
+1 (4) ~ 4) Se a thy (f)— us] ge HE (5 
the magnitude of which 1s determined by the basic motion of the 
missile, i1.e., by the motion of its projection on plane in of the 
nonrotating system of coordinates €ng. The second group 


Bat) 16 (0) — C1 Ret le (2) — aed Ge vane 


eccnsists cf two terms containing, respectively, quantitie. ¢ft: 
and u_(t>. 
& 


Some methcdse of the construction of the current value .f 
function @.ic) abcard the missi.- are examined in the f - 1 owing 


8&6 
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section. As regards the sum (5.3.1), then it, according to its 
structure, is completely similar to equation (2.4.9) for the initial 
ballistic function e(t). Therefore, the approximate construction 
of the sum (5.3.1) can be produced aboard the missile by the same 
methods as were used for the construction of function e(t). Let 

us examine, for example, the construction of this sum aboard the 
missile by means of using readings of two integrators of accelera- 
tions, the axes of sensitivity of which are oriented respectively 
parallel to axes & and n of the nonrotating system of coordinates 
Ent. Similar to equation (2.7.4), the sum (5.3.1) can be represented 
in the form 


Bi = EO— EONS + 9) —w Oe + 
+i — MONS +O — MON + ; 
EO BSE + IN O01 + et) — md 
+18 8S +0) Z. (5.3.3) 


With subsequent conversion of the sum B,(¢) we will use equations 

(2.7.1) and (2.7.2) for tsochronal variations of coordinates and 

projections of the velocity of the missile and also equations of the ; 
expansion in Teylor series (heing limited only to terms of the 

first order of smallness) 


TH—- PE H—Eo) = (0) FO, 


(5.3.4) 
WO— a =n — 16) =) I? , 
U 
and in exactly the same manner b 
“7° a 
é 
ug (0) — ug = uE(t)- us (0% (¢-— 04) yee 
ds: 2, 
eC (5.3.5) 
i 


Taking into nccount, furthermore, that in accordance with 
(2.6.3) 


&7 








du. {a°) opqe 
“$= EF, (5.3.6) 


we transform the expression for sum B,(t) to the form 





Br(t) = a8 (0% Eton + mod + 0 + 
+0—0) [SF 4 SS + 


+R ar (5.3.7) 


Just as in the third chapter, let us disregard in the beginning 
the effect of the change in the acceleration of the force of gravity 
during motion of the rocket not according to the calculated law. 

In accordance with equations (3.2.6), (3.2.7), (3.2.2) and (3.2.3), 
we obtain at small values of the difference t - o* relations 


8E (0) = BSe (0) = Se(t)— SEC) = Se — ISS) + 
FP MON (5.3.8) 
In (t) = BSq(t) = Sa (1)—Sa(l) = S(O — 1S) + 
+(t— a) VoD), 


and also 
“Burg (1) = BV e (1) = Velt) — Ve (t) = Ve) — Velo) + 
+(¢—o')ag(0)), 
bu, (t) = Va - ValN—Vald = Ve(N— IVa (3) + 
+(t—a')a,(6))). (5.3.9) 


Used here again are equations of the exparsion in Taylor series, 

but for projections S; (), S. () of the apparent path of the missile 
and projections re). 1 () of its apperent velocity in the calculated 
mat ton equations (3.2.5) and (3.2.9) are taken into account. 


By means of equalities (5.3.8) and (5.3.9) the sum (5.3.7) 


is reduced to the form 


ur 
ce 


a a 
Put) = S(t) tat Meat Kine — 
[809 He) B+ MCI Z+VG Z]+ 
+ (¢(—3') {(% a _ vees'y] + [SX Sys) __ —Keo s+ 
+ [FP — ey] + [FRE ois] 2 4 8. (5.3.10) 


In the last equation we can produce further simplifications, since, 
according to equality (2.6.1), we have 

FEO — aie = IPO. 1: 

ay . ° e 

Sr a= IE OO: (5.3.21) 


as a consequence, as a result of the integration with respect to 


time within limits of t = 0 to t # o#, allowing for equations 


(3.2.5) 


« 
FE) — vig) = HIE, wd +O, 


Svs) _ Ve) = Valk, (2), VU), ty de pO. (5.3.12) 


Fquations (>».2.i1) and (5.3.12) can be used in the subsequent 
conversion of expression (5.3.10) for the sum B,(t). As a result 
let us obtair. the following final expression for this sum: 


rea 
B= KO) Ft 5, OF 4K, Omg tM oe 
—C+ (6-0) D. (5.3.13) 
Introduced here are designations 
for 


& ‘ 
CH HOV R+ SCOR + VIZ + Oa (5.3.24) | 














Be ete ES a aoa: 
Da BNA WO: al + F ‘ (ECO, msde 


+ AAV OMY: TIE + IED NTT + 
(0), db dy'(0) : 
+[F r+ FAO. CS 9415) 
It is obvious that C and D are constant quantities, which can be 
determined earlier for each specific case of flight of the missile. 


For the construction of expression (5.3.13) aboard the missile, 
one can use the same integrators of accelerations and elements of 
repeated integration as in the solution of the problem on inertial 
control of flight range by the method given in § 2 of the tnird 


chapter. Furthermore, an additional computer and clock mechanism 
are necessary here, 


If, nowever, the control cf flight range of the missile is 
produced with the help of integrators whose axes of sensitivity 
are tarallel to A- and u-directions, then for the fermation of 
expression (5.3.13) there will be required, furthermore, an auxiliary 
device of the type of converter of ccordinates, which continuously 
sclves the system of equations (see § 3 of the third chapter) 


Ve (t)cosd + V,(¢)sin A me Vy {?), 


5 
Vz (t)cosp + V,(é)sinp = V, (¢) (5.3.16) 


relative to quantities ve Ct) ana VG according to data of current 
readings v(t) and vice) above-mentioned intezrators. 


§ 4. Construction of Main Terms of the Function 
of Lateral Deviation 


The second group of corresponding terms of equation (5.2.4) 
for funetion B(t), as was indicated in the previous section, is 
the sun. f two terms 


KO—CIR and we —wy ®. (5.4.1) 
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They can be called the main terms of the function of lateral 
deviation &(t), since they play the main role in the control of 
lateral motion of the missile. ; 


For the construction of current values of the sum (5.4.1) 
aboard the missile, one can use readings V(t) of the integrator 
of accelerations, the axis of sensitivity of which is oriented in 
parallel to axis % of the nonrotating system of coordinates Enz 
(so-called lateral integrator). 


Let us note that, according to equation (1.4.1) (§ 4 of Chapter I), 
there occurs the relation 


we (1) = EO ma) + hes (5.4.2) 


where a, (t) - projection of apparent acceleration of the missile 
on ave &, and fr - projection on the same axis of acceleration of 
the force of the asthe gravity. Thé latter, similar to equations 
(2.5.4), with sufficient accuracy can be represented in the form 


t= — Et. (5.4.3) 


Thus, relation (5.4.2) can be considered the differential equation 
in the desired function ¢(¢t). 


At the instant t = 0 the axes of the starting system of coordi- 
nates zy2, which is connected to the earth, respectively coincide 
with axes of the nonrotating system nt, and the center of masses 
of the missile, on assumption, is found at their common origin. 
Furthermore, at this instant the missile does not have velocity witt. 
respect to the statins system of ecrrdinates ayn. Corseauently, 


initial conditions of equation (5.4.2) are such: 
G0) = 0, BM mu, (0) mus, (5.4.4) 


where ue - velocity of the beginning of the starting system of 
coordinates in the nonrotating system Eng at the instant when 


gL 

















t = 0, i.e., at the beginning of the powered-flight section of 
its flight. 


It is possible to assume that the acceleration of the force 
of gravity proves to have an insignificant effect on the rate of 
change in function C(t). Therefore, as a first approximation it 
is possible not to consider the projection fr in equation (5.4.2) 
and, respectively, consider that 


SEP mac (0- (5.4.5) 


Hence, taking account the second initial condition (5.4.4), we have 


FO ste (t) = Ve (t) + a8. (5.4.6) 
Here 
4 
Ve(td= facnae (5.4.7) 


is the current value of the projection on axis ¢ of the apparent 
velocity of the missile relative to the nonrotating system of 
coordinates Ent. This is the reading of the integrator of accelera- 
ti ns with the axis of sensitivity parallel to the axis %. Secondary 
integration of both parts of equality (5.4.5), allowing alse for 

the first of the initial conditions (5.4.4), leads to the equation 


Ct) = SoC) + uRe, (5.4.8) 


in which 5, (t) - projection on axis ¢ of the vector of the apparent 
path of the missile in the nonrotating systen of coordinates Ent. 


In accordance with the last equation, for the construction 
aboaurd the missile of current values ((t), an additional integration 
-f peadings of the aforementioned integrator of accelerations 1s 


i 


necessary. Actually 


t 
Set) = Ve (nat (5.4.9) 
® 
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On the basis of equations (5.4.6) and (5.4.8) the expression 
for the sum (5.3.2) is reduced to the form 


‘ 7 2, ian “ «, 2 
Ba (t) = 15e (8) Fel — OV Ge EMCO Fe Ed Gar (5.4.10) 


This sum, just as sum B,(t)s obtained in § 3 of this chapter, 
can be constructed directly aboard the missile. 


Thus, under the assumptions made about the unimportance of 
quantity fy it is possible by means of integrators of accelerations 
and computers to obtain current values of the function of lateral 
deviation ° 


Be) = Bit) + Ba(!). (5.4.11) 


§ 5. Calculation of the Acceleration of the Force 


of Gravity in the Construction of the Function 
of Lateral Deviation 


Methods given in the fourth chapter, allow, as will te shown 
below, approximately integrating equation (5.4.2) for function 
c(t) taking into account the projection of the acceleration of the 
force of the earth's gravity fre 


Let us note, first of all, that with an accuracy sufficient for 
practice at not too great an extension of the powered-flight section 
of the flight of the missile in equation (5.4.3), the distance po fror 
the missile to the center of the earth can be substituted by the 
radius of the earth R. After this let us arrive at equation 


f= — v0, (5.5.1) 


in which quantity v (Schuler frequency) is determined by equation 
(4.3.4). Finally equation (5.4.2) can be replaced by the following: 


FE on ag (0) — v (4), (5.5.2) 
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Let us integrate witis respect to time the right and left sides of 
this equation within limits of 0 to ¢ taking into account the second 
initial condition (5.4.4). Equating the obtained results, we arrive 
at integral differential relation 


: . 
Pere tena + ue. (5.5.3) 


Let us prceduce the same operation of integration agair but taking 
inte account the first initial condition (5.4.4). As a result we 
ottain the integral relation 


Be KO I Po un|ee + wpe, (8.5.4) ‘ i 


whicn is the equation for the determination sf lateral shift ef the 
missile ¢ft) according to the assigned function S,(t). The expresslon 
represented by equation (5.4.8) can be considered the first approxi- 
mation to the solution of this equation. To obtain the following 
acproximatiorn, in which already tne effect of acceleration of the 

force of gravity will be taken into account, it is possible to use 

the method given in the fourth chapter. First of all, let us assume 

in relation (5.5.4) t # 3. Then it takes the form : 


2 ter nes moaning ctetnd 


et 
(2) = Secay— vt (nae + up. (5.5.5) 


Let us further replace in the right side of equation (5.5.4) function 


a(t) by its approximation in the form of a polynomial of the second 
power 


(2) — & 
Gey — une + COLES os, (5.5.6) 
the coefficients of which are selected so that initial conditions 


(5.4.4) would be satisfied, and, furthermore, in order that the ‘ 
polynomial would take value ¢(o) at the instant ¢ * o of termination : 


yd 


& 
= 
z 
* 





3 
3 
3 
4 
3 
3 
3 
Z 
3 
5 
g 
4 
3 
3 
3 
3 
j 
g 





of the powered-flight section of real m-ticn of the mizssile.! 


if then we produce a similar replacement alsc in relation 
(5.5.3) and accomplish the necessary operations of integration, 
then we will arrive at equality 


uz (3) = Ve (9)— viug v% (0) 5 + ue, 


a a (5.5.7) 
" 6(0) = Sz (5) — v8ue > — vb (3) a5 + U8s, 


which can be examined as the equation for determining the value of 


the function itself Z(a) and its time derivative u,(o) at the 
instant t = o. 


The solution of the second equation (5.5.7) can be obtained 
with sufficient accuracy if in its right side quantity ¢(o) is 
substituted oy its first approximation, according to equation 
(5.4.8), i.e., if we assume 


(6) = Sz (5) + ups. (5.5.8) 
as a result we obtain equation 
£(0) = (1— SF) 5, (a) + (1 — Ft) ous. (S.8.0° 
tr 


We produce the same replacement in the first equation of (5.5.7). 
t 


Finally we arrive at the following approximate equation for the 


_derermination of the project#en of the actual velocity of the 


missile on axir 3 


tg (9) = Ve (3) — Se (3) + (1 — 7) (5.5.20) 
Considering the argument o to be close to of, it iz ,orsible 


Withent great loss of accuracy in equations (5.5.9) and ‘°.5.10) 


bovoare ts a possibility, similar to § 1 of the fourth shapter, of 
an approximation by means of a polynomial of the third nower fr 
which the derivative when t = G still turns into 4, (o}. 








to substitute argument o in terms containing the small factor se, 
by the rated value of duration of the powered-flight section o*. 
Substituting in other places of these equations argument o for the 
current time t, we obtain the following final expressions for c(t) 


and u,(t) 


t= [4 aes eer] s (+ (1 se wer] tut, 


te (1) = Ve (2) — 7 Sec) + [1— SSF we, (5.5.11) 


Equations (5.5.11) lead to a presentation of the sum Batt), more 
accurate in comparison with expression (5.4.10), namely: 


a(t) = ft — “SF sect — 2 Jan — By 


+{hO- Es. +[ _ Oe] up ug} 2. oibere7 


The construction of function Bo(t) aboard the missile can be 
preduced by means of the integrator of acceleration with axis of 
sensitivity parallel to axis ¢ (and also, of course, the integrating 


device). 


in this bock, as was already Indicated in the intreducticn, an 
necount of problems or the theerv of inertial sortro? o8 the flight 
of ballistic ini:siles was produced under the i:sunption of the ideal 
stabilization of axes of sensitivity of ‘awtc-neters and integrators 
of accelerations and completely accurate cper:tton of the latter, 
i.e., accurate Indication by them of current values of projections 


of the veetar of apnarent acceleration of the missile and interrals 


of th ce Vallure, 


An estimate of the effect of instrumental errors of gyrcscopic 
deviees ami dnverrators of accelerations on the marnicude of possible 
dewlatiar (os the adestle with Plight at the ussimicd ts. st and also 





an analysis of causes of the appearance of these errors and develop- 
ment of measures on their elimination represent an independent 
interest and fall outside the framework of this book. 


Let us note, finally, that problems of inertial control of 
space missiles can be solved by the same methods as the problem 
examined above about the flight of ballistic missiles within limits 


of the earth. 
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APPENDIX 


DERIVATION OF EQUATIONS FOR COSINES OF ANGLES BETWEEN ; 
AXES OF TWO SYSTEMS OF COORDINATES TURNED RELATIVE : , 
TO ONE ANOTHER AT A FINITE ANGLE 


In manuals on analytical geometry and according to theoretical 
mechanics, equations for the cosines of angles hetween axes of two 


systems of ccordinates arbitrarily turned relative to each other are 
derived usually with the help of several additional geometric con- 
structions. The derivation of these equations is also possible 

cy means of the frequent use of tables of cosines of angles between 
axes of the basic and auxiliary system of coordinates, each of 
which is turned relative to the previous around one of the coordi- 
nate axes’ common with it. Specifically, the main axis of rotation 
is simultaneously the axis of coordinates of two auxiliary systems 
turned one reiative to the other at the same angle as the main 


cee OI BR aemeciamannmee ne HO 





axes. 


Given below is a purely analytical derivation of the mentioned 
equations, founded exclusively on the simplest theorems of 
analytical geometry, and some of these theorenm- are used in 


vectorial form. { 


let us assume that the system of coordinates €nt 1s turned at 
an arbitrary angle $ around a certain axis d, which passes through 


‘See, for example, Whittaker. Analytical Kinetics. M.-L., 1937 
and A. Yu. ishlinskty. Mechanics of Gyroscopic Systems. IJzd-vo 
AN SSSR, 1963. 
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the origin of this system.’ The final ;.osition of the system is 
designated by xyz. If 1, m and n are -osines of angles which form 
the mentioned axis of the final turn d with axes of the system of 
coordinates €nt, then it is obvious that the same magnitudes are 
respectively cosines of angles between the axis d and axes of the 


system xyz. Specifically, we have 
i = cosa, (1) 


where a@ - angle which forms the axis of rotation d with the axis € 
of the system of coordinates €ng and simultaneously the same axis 
of rotation d with the x axis of the system zyz (see Fig. 14). 


Fig. La. 





Let us designate the direction cosines of the # axis relative 


te the system of coordinates Eng respectively by a, 2 anc e. They 


are in this case the desired quantities. Besides the apparent 


equality 
a+ oipctat — (2) 


‘In §§ 2 and 3 of the second chapter of this book such an axis 
vais the axis of retation of tie earth, and angle ¢ was the product 
of the angular velocity of the earth U for time t, which passes fror 
that instant when the axes of systems of coordinater Enz and xyz, 


rospeetively, coincided. 


go 





the direction cosines ae, b and e sacisfy also the relation 


la +mb+ne= I, (3) 


which, in accordance with the known equation of the analytical 
geometry, expresses the magnitude of the cosine of angle a between 
the x axis and axis of rotation d. 


For the determination of the three desired quantities a, b and ec, 
one equation is necessary. Specifically, it can express the fact 
that plane xd is turned relative to plane Ed about the axis of 
rotation d at the assigned angle ¢. Thus the measure of the dihedral 
angle between the mentioned planes also should be equivalent to 


angle ¢. 


Sines the dihedral angle is measured by the angle between two 
rerpendiculars to the planes forming it, then we will introduce 
two vectors h and k, the first of which is perpendicular to the 
niane Ed, and the other - plane zd. According to the property of 


the scalar pr. :vet of two vectors, we have 


hk cos @ = heke + hake + hike, (4) 
where hes has he - projections on the axis —, n, % of vector h 
and, respectively, Kes ky ky ~ projections on these axes of vector 


k, 


Let us take as vector hk che vectoria: przduct of the single 
vector e. located on the € axis, and singl* vector ges which 
possesses the direction of the axis of rotation d. Projections of 
the first of them on axes €, n and ¢, of course, are equal, 


respectively, to numbers 


1, 0, 0, (5) 
amd of the second 

i, m,n. if) 
According fo ocato of the comperitlon of the vectorial: rreduct, 


we obbala. 


o~ om og 
3 om 


& 
0 |: (7) 
a 








and, consequently, 


dg =0, Aym—nr, Ag = m. (8) 
In turn for vector k we take the vectorial product 


taaxoe[e (9) 
i 


3 4 


f 





Here the second line of the determinant is made up of projections 
of the single vector 2° on the axes €, n and ©, equal, of course, 
to the direction cosines of the x axis in the system of coordinates 


Eng. 


Expanding the mentioned determinant, we arrive at the following 
equality: 


byand—me, Rye le—na, hy = ma— Ib. (10) 


Moduli of both vectorial products (7) and (9) are identical 
and equal to the sine of angle*a between axes — and d or, which is 
the same, between axes x and d. Thus, taking into account still 


another equation (1), we have 


kmh wsina = VI —P. (11) 


By means of equalities (8), (10) and (11) relation (4) can be 
presented in the form 


(1 — &@)cosg = — I (md + ne) 4 (a? + m)a. (12) 


It is the deficient third equation, besides equations (2) and (3), 
for the search of the three unknown cosines a, b and c, which 
determine the direction of the x axiz with respect to axes of the 


avstem of coordinates €nzt . 
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According to equation (3), we have 
mb + new if —a). (13) 

Furthermore, 
m*+at ai — FB. (14) 


Taking account the last two equalities in equation (12), we obtain 


after reducing similar terms to the equation for one of the desired 
quantities 


eo = cose + (i — cos 9). (15) 


lo determine quantity b, we exclude from relation (2) cosines 


a and ¢ by means of equations (12) and (15). As a result we arrive 
at the quadratic equation 


IP + (1—P)cose] +4 SP + 
+ (1 — 7) cos 9] 1 — bua}® = f, (16) 


whicn after simplifications with the help of relation (14) is 
reduced to the form 


5-25 im (8 cos q) +Pm* (1—cos @)*—n* sin o=0. (17) 

Of the two roots ef this equation 
4, = bn (1 — coag) +a clog, 3 
& = im (i — cos @) — naing or) 


.e snould diiecuss the first. Actually, assuming specifically, 


lem=0, axl, =F, (19) 
we cucain accordingly 
k= +1, b= — 1. (20) 
10e 
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However, taking into account the accepted designations (see Fig. 


14), 
we have in this case 


5 = cos rn = +1 (21) 


and, consequently, the second root of equations (18) should be 
dropped. Thus, in general 


Substituting now expressions (15) and (22) for @ and b into equation 


(3), we will arrive at the following equation for the determination 
of quantity e, namely: 


e = ln (1 — cosg) — m sing. (23) 


Similarly cosines of angles between the y axis and axes &, 
n, 6 and further between the 2 axis and the same axes €, n, & are 
dgetermined. Finally, the table of cosines of the angles between 


axes of the system of coordinates xyz and Eng can be represented 
in the following form: 


t " C 
#(i—cose@)M+4 ({—cos¢)ml4+ (1—cos@)ni— 


+ co3¢ + nsing —msing 
y (1 — coe ¢) Im— {t—cos9)m*+ (1—cosq) mn + 
—ansing + cos@ +ésing 
3 (1—cos@)in+ (1—cos@) mn— (1 — cos@) nt + 
+msing —tlsin¢g +cos@ 24) 


"his table was used in the second chapter 2f this t x0k. 
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